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13  abstract  iMai.T 200 

V  ^ 

During  the  past  year  work  at  Los  Aleimos  concerning  early  time  nuclear  phenomenology  at 
high  altitudes  has  involved  a  number  of  diverse  problems.  These  include:  (1)  simulations 
of  magnetic  field  and  debris  dynamics  in  VHANE  bursts,  (2)  calculations  of  the  energy 
deposition  spectrum  in  HANE  bursts,  (3)  study  of  debris  characteristics  in  both  VHANE 
and  HANE  regimes,  (4)  development  of  a  simple  debris  slowing  and  deposition  model 
leased  on  electromagnetic  instabilities  due  to  debris  streaming  along  the  magnetic  field, 
and  (5)  some  preliminary  calculations  of  multiburst  effects.  ,The  principail  results  are 
first  summarized,  and  then  the  details  of  the  individual  caJcuIptions  are  discussed  in  the 
following  five  separate  sections. 
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Summary 


We  briefly  summarize  work  carried  out  in  FY91  at  Los  Alamos  in  relation  to  high  altitude 
nucleEir  bursts.  This  summary  is  subdivided  into  five  parts,  each  representing  different 
aspects  of  the  work.  Details  of  the  individual  tasks  are  fovmd  in  the  following  sections  of 
the  report,  each  of  which  is  self  contained.  For  each  part  we  briefly  discuss  the  purpose 
of  the  work,  how  it  was  carried  out,  the  principal  results,  the  consequences,  and  proposed 
future  studies.  The  sections  are  titled:  (1)  VHANE  simulations,  (2)  Energy  deposition 
spectrum  in  a  HANE  burst,  (3)  Debris/air  interactions,  (4)  Parallel  streaming  deposition 
model,  and  (5)  Multibursts. 

1.  VHANE  simulations.  This  study  involved  determining  the  size  of  the  magnetic  cavity 
in  a  VHANE  burst  as  well  as  characterizing  the  properties  of  the  debris  and  air  ions  after 
the  magnetic  field  recovered.  We  used  two  hybrid  (particle  ions,  massless  fluid  electrons) 
simulation  models  independently  to  investigate  how  the  plasma  expands  across  the  mag¬ 
netic  field.  A  newly  developed  hybrid  code  that  includes  electromagnetic  radiation  wels 
employed  to  study  2-D  (r-z)  debris  expansions.  It  was  fovmd  that  the  cavity  extended 
out  to  about  Rb  —  O.8R3,  where  R3  is  the  classic  magnetic  confinement  radius  obtained 
by  equating  the  initial  debris  kinetic  energy  to  the  magnetic  field  energy  in  a  sphere  of 
radius  R3.  It  was  cdso  shown  that  the  cavity  size  was  independent  of  the  background  air 
density  (riair  ‘C  ndebri$  or  riair  =  0).  A  standard  hybrid  code  was  also  used  to  study 
VHANE  bursts  expanding  in  the  two  directions  perpendicular  to  the  magnetic  field  and 
to  characterize  the  debris  when  the  magnetic  cavity  collapsed.  The  final  debris  radius 
was  about  0.5iZ2,  where  Ri  is  the  classic  confinement  radius  for  a  2-D  expansion.  Some 
ions  remained  a  distsmce  1  —  2p,  {pi  =ion  gyrorEwiius)  outside  of  the  main  debris  cloud; 
this  fraction  increased  with  pi/Ra-  The  debris  remained  hot,  with  a  final  temperature 
Td  O.5E0  {Eo  =initiaJ  debris  kinetic  energy),  with  most  of  the  thermalization  occurring 
as  the  magnetic  field  recovered.  The  debris  was  also  structured,  with  flute  modes  forming 
on  the  surface  at  the  end  of  the  expansion,  as  seen  in  numerous  other  calculations.  The 
structure  persisted,  as  the  debris  ions  recollapsed  with  the  magnetic  field.  While  the  cal¬ 
culations  were  done  in  2-D,  one  anticipates  similar  results  in  3-D,  but  with  a  smadler  final 
radius  (due  to  elongation  of  the  plasma  along  B)  and  with  the  debris  ions  about  a  third 
cooler  (energy  partitioned  into  three  rather  than  two  degrees  of  freedom).  Overall,  these 
results,  along  with  those  obtaiined  by  MRC,  BRA,  and  LLNL,  have  led  to  a  very  good 
understanding  of  the  magnetic  field  and  debris  properties  in  a  VHANE  burst.  Except  for 
multiburst  effects  (see  Part  5),  no  new  effort  in  this  area  is  needed  or  plannned. 

2.  Energy  deposition  spectrum  in  a  HANE  bvirst.  This  study  was  not  part  of  our  planned 
program  for  the  year,  but  was  carried  out  to  help  resolve  a  controversy  between  the  BRA 
and  MRC  deposition  spectra.  Thompson’s  use  of  these  spectra  in  NORSE  calculations 
yielded  significant  differences  in  the  phenomenology,  leading  to  a  discussion  of  the  ori¬ 
gin  of  these  spectra.  We  carried  out  some  2-D  hybrid  simulations  of  relatively  low  yield 
HANE  bursts  (i.c.,  initial  debris  Mach  number  of  16)  to  obtain  spectra  of  energized  air 
ions.  As  these  spectra  were  calculated  locally,  rather  than  at  lower  altitudes  where  the 
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deposition  actually  occurs,  they  were  not  deposition  spectra  per  se,  but  exhibited  char¬ 
acteristic  features  of  kinetic  ion  {e.g.,  BRA)  calculations.  The  features  included  a  double 
peaked  spectrum  in  V||  and  a  maximum  Uji  >  Vb  (=initial  debris  expansion  velocity).  We 
2dso  showed  examples  of  spectra  from  coUisionless  shocks  in  space  that  exhibited  similar 
features.  While  supporting  the  basic  BRA-type  deposition  spectrum,  this  work  has  also 
suggested  the  need  for  further  BRA/MRC  spectrum  comparisons  in  NORSE  calculations 
(J.  Thompson  has  zdready  carried  out  some  additional  nms)  to  further  vmderstand  what 
aspects  of  the  phenomenology  are  sensitive  to  what  features  of  the  spectrum.  As  the  ki¬ 
netic  simulations  also  tend  to  show  that  a  significant  fraction  of  the  debris  energy  remains 
undeposited,  there  is  also  a  need  to  imderstsmd  where  this  energy  eventually  goes  and  how 
to  model  it  appropriately  in  the  systems  codes.  Finally,  the  question  of  the  importance  of 
the  deposition  spectrum  suggests  the  use  of  the  AGNES  laser  facility  to  'lirpctly  measure 
debris  energy  spectra  (both  as  a  function  of  space  and  time)  and  compare  with  appropriate 
numericed  calculations. 

3.  Debris/air  interactions.  This  work  involved  a  series  of  hybrid  simulations  in  the  HANE 
and  VHANE  regimes  to  characterize  the  properties  of  the  debris  and  air  ions  as  a  function 
of  debris  Mach  number.  Mo-  The  calculations  were  done  in  2-D,  perpendicular  to  the 
magnetic  field.  Generally  it  was  found  at  low  Mo  (VHANE  regime)  that  the  debris 
remained  well  confined  and  heated  as  the  magnetic  field  collapsed.  A  small  number  of  air 
ions  were  trapped  in  the  debris  during  the  collapse  of  the  field  and  were  likewise  heated. 
At  higher  Mo  (^4),  a  collisionless  shock  formed  that  heated  the  air  ions.  The  debris  did 
not  completely  collapse  in  this  case,  because  the  present  simulations  did  not  realistically 
let  the  magnetic  field  recover.  At  Mo  =  aspects  of  both  low  and  high  Mo  expansions 
were  observed  to  occur:  structuring  of  the  debris  and  resultant  heating  as  the  magnetic 
cavity  collapsed,  like  the  low  Mo  expansions,  and  the  formation  of  a  weak  shock  and 
heating  of  the  air  ions,  as  in  high  Mo  expansions.  There  was  no  noticeable  large  coupling 
to  magnetosonic  waves  at  Mo  =  1,  as  some  EMP  models  predict.  As  a  function  of  Mo, 
the  final  debris  radius  increased,  while  the  energy  in  the  confined  debris  and  air  decreeised, 
cis  Mo  increased.  These  calculations  have  shown  that  it  is  possible  to  characterize  the 
properties  of  the  debris  and  air  ions  when  the  debris  has  stopped  and  the  magnetic  field 
recovers.  At  higher  Mach  number,  some  improvements  in  the  calculations  au’e  needed  to 
follow  the  evolution  further  in  time.  However,  the  preUminary  results  here  axe  encouraging 
and  provide  the  basis  for  developing  a  complete  model  for  the  debris  location  zind  energy 
spectrum  in  HANE  and  VHANE  bursts. 

4.  Parallel  streaming  deposition  model.  This  study  involved  the  development  of  a  model 
for  the  slowing  and  deposition  of  debris  ions  streaming  relative  to  the  air  along  the  magnetic 
field  via  electromagnetic  ion  beam  instabilities.  We  reviewed  properties  of  the  instabilities 
and  developed  simple  criteria  for  their  application  to  VHANE/HANE  bursts.  To  apply 
these  criteria  we  considered  two  simple  models  for  the  debris  expansion.  The  first  model  as¬ 
sumed  that  the  debris  expanded  uniformly.  In  this  case  the  debris  density  quickly  dropped 
below  that  of  the  air,  and  the  resonant  ion  beam  instability  provided  the  coupling.  This 
model  indicated  that  the  instability  was  reduced  at  lower  altitudes,  as  the  debris  density 
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faJls  off  rapidly  from  the  biirst  point,  so  any  coupling  that  might  occur  happens  at  higher 
altitude.  The  second  model  assumed  that  the  debris  remained  in  a  narrow  coupling  shell, 
and  hence  the  debris  density  exceeded  that  of  the  air  down  to  low  (~  200  km)  altitudes. 
In  this  case  the  nonresonant  ion  beam  instability  was  the  dominant  coupling  mechanism, 
and  the  stopping  distance  due  to  anomalous  wave-particle  collisions  was  shorter  by  about 
three  orders  of  magnitude  than  that  due  to  binary  collisions.  The  two  models  in  some  sense 
represented  extreme  conditions,  with  the  actual  situation  somewhere  in  between.  Never¬ 
theless,  the  cedculations  do  suggest  that  anomalous  debris-air  coupling  by  electromagnetic 
ion  streaming  instabilities  can  lead  to  the  deposition  of  the  debris  ions  at  somewhat  higher 
altitudes  than  is  usually  assumed.  Thus  far,  the  calculations  have  been  limited  to  point 
by  point  comparisons  of  anomalous  versus  binary  coUisional  effects.  What  is  needed  is  a 
model  for  the  continuous  slowing  of  the  debris  ions  by  anomalous  effects,  allowing  for  the 
spread  of  the  debris  and  the  increase  in  the  air  density  at  lower  zdtitudes.  Also  needed  are 
simulations  of  such  spatially  varying  situations  to  verify  the  assumptions  of  the  model. 

5.  Multibursts.  The  final  topic  we  addressed  was  the  question  of  new  effects  due  to  multiple 
bursts.  We  carried  out  both  full  particle  and  hybrid  simulations  of  two  simultaneous  bursts 
in  the  VHANE  regime.  We  chose  the  VHANE  regime,  because  most  of  the  physics  of 
a  single  bvurst  is  now  well  understood.  We  foimd  that  many  aspects  of  the  multiburst 
calculations  were  similar  to  the  single  burst  case.  We  observed  the  initial  expansion  of  the 
debris  ions  and  thier  compression  into  a  thin  shell,  which  ultimately  structmed.  However, 
the  two  simultaneous  bursts  gave  rise  to  greater  heating  of  both  the  ions  and  electrons 
than  in  a  single  burst,  and  structure  tended  to  be  magnified  rather  than  suppressed.  We 
interpreted  these  results  as  being  due  to  the  interaction  of  the  plasma  from  one  burst  with 
the  macroscopic  and  microscopic  fields  of  the  second  burst,  rather  than  due  to  some  sort 
of  ion /ion  instability.  Such  instabilities  would  tend  to  be  numericzdly  suppressed  in  the 
present  calculations.  These  preliminary  results  suggest  that  treating  multiple  bursts  as 
superpositions  of  single  bursts  is  likely  to  underestimate  heating  and  structuring  effects, 
perhaps  by  a  large  amount.  Even  two  bursts  will  lead  to  a  larger  diamagnetic  cavity, 
enhanced  fluctuations,  and  increased  (factor  of  two)  heating.  It  is  felt  that  this  is  an  area 
where  the  AGNES  facility  coiild  be  most  useful.  Experiments  involving  two  explosions 
could  be  carried  out  in  a  straightforward  manner  to  study  these  types  of  enhancements  as 
well  as  to  investigate  ion/ion  coupling  in  a  coUisionless  medium. 
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1.  VHANE  Simulations 


1.1  Introduction 

Simulations  of  VHANE  bursts  have  been  carried  out  with  two  hybrid  algorithms,  one 
being  a  “standard”  hybrid  model  that  requires  a  background  plasma  and  mzdces  the  low 
frequency  (Darwin)  approximation  and  the  other  being  a  new  model  that  allows  expansion 
into  vacuum  and  keeps  ah  electromagnetic  radiation  effects.  The  cadculations  have  been 
done  to  determine  the  maximum  spatial  extent  of  the  debris  as  weU  as  the  debris  location 
and  energy  distribution  after  the  magnetic  bubble  coUapses.  We  first  describe  the  results 
with  the  new  hybrid  algorithm  in  Section  1.2  and  then  calculations  with  the  standard 
algorithm  concerning  the  debris  characteristics  in  Section  1.3.  Conclusions  are  given  in 
Section  1.4. 


1.2  Debris  Extent  Calculations  with  a  New  Hybrid  Model 


1.2. A  Background 


In  the  VHANE  regime  the  expanding  debris  is  contained  by  the  Earth’s  magnetic  field.  The 
magnetic  confinement  radius  Rz  is  easily  estimated  by  equating  the  initial  kinetic  energy 
of  the  debris  (ion  mziss  mo,  number  Nd,  velocity  Vb)  with  the  equivalent  magnetic  field 
energy  in  a  volume  of  radius  i?3, 


to  give 


1  .^2  ^TT  3 

^NomoVi,  =  —-R, 


R,  =  (?.ypyy>P)./3 

D- 


(1.1) 


(1.2) 


If  one  instead  assumes  a  2-D  expansion  perpendicular  to  B  (as  we  will  later),  one  finds  a 
similar  formula, 

fi,  =  (1,3) 

with  Nd  the  debris  line  density  (particles  per  cm^). 


The  key  issue  is  whether  in  an  actual  3-D  expansion,  the  debris  extends  out  to  R3.  Kilb 
[1990],  for  example,  argues  that  one  should  use  1.5B  rather  than  B  in  Eq.  (1.1)  so  that 
the  actual  debris  expzmsion  radius,  here  called  Rb,  reduces  to 


Rb 


ANottidVI 
^  3B2  ’ 


0.76^3 


(1.4) 


Alternatively,  Gisler  and  Lemons  [1989]  have  concluded  that  because  the  debris  is  not 
expanding  everywhere  perpendicular  to  the  magnetic  field,  as  assumed  in  Eq.  (1.1),  Rb  is 
smaller  by  a  factor  of  (0.5)^/^  so  that 


Rb  ~  0.88723 


(1.5) 
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For  typical  VHANE  conditions  (iVomo  ~  10® gf,  B  ~  0.3G,  Vb  ~  2000km/sec),  one  finds 
R3  ~  llOOkm.  Thus,  the  difference  between  Eqs.  (1.1)  and  (1.4)  or  (1.5)  for  the  maximum 
extent  of  the  debris  is  about  260  km. 

1.2.B  New  Hybrid  Algorithm 

Particle  simulations  that  we  have  carried  out  [Winske,  1990],  as  well  as  hybrid  simulations 
by  Brecht  [private  commvmication],  Simonson  and  Hewett  [1991],  and  us  generally  tend  to 
show  Rg  ~  O.SiZa,  in  agreement  with  Eq.  (1.4),  except  in  a  few  cases.  To  help  resolve  the 
differences,  we  have  used  a  new  hybrid  algorithm  [Jones  et  ai,  1989],  which  can  treat  debris 
expansion  into  a  magnetized  vacuum  as  well  as  into  a  low  density  background  plasma  and 
does  not  make  the  usual  low  frequency  approximation  (hence  includes  electromagnetic 
radiation). 

In  a  “standard”  hybrid  code,  one  assumes  massless  electrons  (m*  =  0)  and  solves  the 
electron  momentum  equation  for  the  electric  field 

nienedVe/dt  =  0  =  -ent(E  +  K  x  B/c)  —  Vp*  +  encTjf  (1.6) 

(where  is  the  electron  fluid  velocity,  E  is  the  electric  field,  B  is  the  magnetic  field,  pt 
is  the  electron  pressure,  rj  is  the  resistivity,  and  J  is  the  current).  The  electron  current 
(Jg  =  —eueVe)  is  obtained  from  the  ion  current  and  Ampere’s  law  (in  the  Darwin  limit) 

Jt  —  ~Ji  +  7"^  X  B  (1.7) 

47r 


The  new  algorithm  solves  MeixwcU’s  equations  as 


=  cV  X  B  —  4n{Ji  +  Jj) 


dB  -  - 

^  =  -cVxE  (1.9) 

In  this  case  Eq.  (1.6)  is  solved  for  and  the  Eqs.  (1.8-1. 9)  are  advanced  as  in  a  standard 
explicit  electromagnetic  pzu-ticle  code  [Birdsall  and  Langdon,  1985]. 

1.2.C  Simulation  Results:  No  Backgroxmd  Plasma 

In  order  to  show  the  utility  of  this  approach,  we  have  considered  the  axisymmetric  (r-z) 
expansion  of  a  debris  cloud  into  a  vacuum,  in  which  is  embedded  a  uniform  magnetic  field, 
B  =  BoZ.  Instead  of  expressing  the  initial  parameters  in  terms  of  code  units,  we  can 
express  lengths  in  terms  of  B3  and  the  time  in  terms  of  r  =  Vbt/Ba.  The  initial  debris 
cloud  has  a  cosine-shaped  density  profile  that  extends  out  to  0.125  B3.  The  ions  are  cold 
with  directed  radial  velocity  Vg  such  that  the  ratio  of  the  gyroradius  {pi  =  VofCli,  where 
n,  =  eBo/niiC  is  the  ion  gyrofrequency)  to  R3  is  0.0625. 
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Figure  1.1  shows  a  time  history  of  the  debris  ion  kinetic  energy  Ed  and  the  magnetic  field 
energy  Eb  (in  terms  of  the  initial  debris  energy).  As  expected,  Ed  decreases,  while  Eb 
increaises  accordingly.  Eventually,  Ed  reaches  a  minimum  as  the  debris  stops,  while  Eb 
maximizes  as  the  cavity  reaches  its  largest  radius;  this  occurs  at  r  ~  1.  Figme  1.2  shows 
a  radial  cut  of  the  magnetic  field  at  r  ~  0.9.  The  field  is  completely  excluded  from  the 
cavity,  which  extends  out  to  Rb  ^  O.TSRs.  At  the  edge  of  the  cavity  the  magnetic  field 
is  compressed  to  a  value  about  40%  above  the  ambient  field.  The  magnetic  field  diffuses 
to  some  extent  (the  applied  resistivity  in  terms  of  a  collision  frequency  is  10) 

so  that  the  magnetic  field  (and  density)  scale  length  is  about  pi.  Figure  1.3  shows  a 
contour  plot  of  the  magnetic  field  at  r  ~  0.9.  One  sees  a  smooth  cavity,  with  an  elongation 
factor  (Lz/Rb  ~  1.4),  consistent  with  theoretical  expectations  [Gisler  and  Lemons,  1989]. 
Finally,  Figure  1.4  shows  the  debris  ions  at  various  times.  The  ions  compress  into  a  shell 
as  they  expand,  eventually  reaching  Rb  —  O.lbRi.  A  small  number  of  ions  zire  beginning 
to  leak  out  the  ends  of  the  cloud  along  the  z  axis  at  this  time. 

1.2. D  Calculation  with  Background  Plasma 

Brecht  [private  communication]  had  suggested  that  differences  in  these  new  hybrid  calcu¬ 
lations  and  in  our  previous  particle  simulations  with  some  of  the  hybrid  calculations  done 
at  BRA  and  LLNL  might  be  due  to  the  absence  of  a  background  plasma  in  our  work.  The 
above  calculation  was  thus  repeated  including  a  low  density  background  plasma.  In  the 
previous  calculations  the  initial  debris  density  was  n^/nc  =  17280  {rif.  =code  reference 
density),  with  an  expansion  velocity  of  Vd(va  =  0.17,  with  va  =Alfven  speed  in  terms  of 
the  reference  density.  Here,  the  backgroimd  density  was  no/”c  =  3-75,  so  that  the  debris 
would  overrim  a  backgroimd  mass  that  was  10%  of  the  debris  mass  if  it  expanded  to  R3. 
Figure  1.5  shows  a  contour  plot  of  the  debris  density  at  r  ~  1.0.  Again,  the  ions  extend 
out  to  Rb  —  0.75/23,  consistent  with  the  earlier  result.  Also,  no  strong  anomeilous  heating 
due  to  the  interaction  of  the  debris  ions  with  the  background  plasma  was  observed. 

1.3  VHANE  Debris  Characteristics 

We  next  consider  two  dimensional  expansions  in  the  plane  perpendicular  to  the  magnetic 
field  to  determine  properties  of  the  debris  after  the  bubble  collapses.  The  2-D  calculations 
allow  the  development  of  flute  intabilities  on  the  surface  of  the  cloud  and  provide  a  good 
comparison  with  the  r-z  simulations  of  this  problem  by  Brecht  [private  communication] 
and  Simonson  and  Hewett  [1991]. 

1.3.  A  Description  of  One  Simulation 

In  this  case  we  use  a  standard  hybrid  algorithm  [e.g.,  Winske  and  Quest,  1986],  making 
the  usuzd  assximptions  of  low  frequency  aind  miissless  electrons.  The  debris  has  an  initial 
Gaussian  shaped  profile  (n£j(r)  =  nDexp(--r^/r^),  =  O.I6/22)  with  expansion  velocity 

Vd  '>  the  ion  thermzd  speed  (uj,  expressed  in  terms  of  =  SnrioTDlBl,  Td  =  0.5m£)i;?) 
The  initial  peak  debris  density  is  riDlno  =  41  with  mDlrrio  =  4,  and  =  0.01.  The 
expansion  velocity  expressed  in  terms  of  a  debris  Alfven  Mach  number,  is  Md  =  Vd/va  = 


6 


B./B 


Figure  1.2.  Radial  cut  of  the  magnetic  field  at  r  'w  0.9,  showing  the  cavity,  the  B  field 
enhancement  on  the  outside,  and  the  steep  gradient  at  the  edge  of  the  cavity. 
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Figure  1.3.  Contour  plot  of  the  magnetic  field  at  r  '>«  0.9,  showing  smooth  cavity  at 
maximum  expansion,  somewhat  elongated  along  B. 
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Figvire  1.5.  Contour  plot  of  debris  density  at  r  =  1,  showing  expansion  to  same  radius  in 
a  calculation  with  a  weak  background  plasma. 
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0.5,  with  va  defined  in  terms  of  the  background  density  and  ion  mass.  Figure  1.6  shows 
the  debris  density  at  various  times.  The  debris  expands  out  to  a  radius  r  ~  722-  Flute 
modes  begin  to  develop  as  the  debris  stops.  The  fiutes  elongate  as  the  plasma  collapses, 
leading  to  a  rather  well  contained,  highly  structured  debris  cloud.  Figure  1.7  displays  3-D 
perspective  plots  of  the  debris  density,  the  debris  energy  {Ed  =  0.5m d[vI  +  u^]),  and  the 
background  air  ion  energy,  when  the  cloud  has  collapsed  =  80).  The  highly  structured 
nature  of  all  three  quantities  is  evident.  Figure  1.8  shows  a  radial  slice  of  the  debris  density 
and  energy.  We  characterize  the  collapsed  cloud  in  terms  of  an  inner  radius  rj  and  an  outer 
radius  r2 ,  along  with  corresponding  average  values  of  the  debris  density  and  debris  and  adr 
energies  in  the  regions  r  <  ri  and  rj  <  r  <  ra.  ri  is  defined  as  the  radius  where  the  debris 
density  drops  to  20%  of  its  peak  value,  while  ra  is  the  radius  of  maximum  extent  of  the 
debris.  This  is  not  the  only  characterization  possible,  but  does  allow  one  to  then  compare 
these  properties  as  a  fimction  of  burst  parameters. 

In  Figure  1.9  we  plot  the  energy  distribution  f{E)  (E  =  0.5m£)[uj -t-u^j),  normalized  to  the 
initial  debris  energy  Eo,  for  the  debris  ions  at  Hit  =  40  [top  panel]  at  the  time  of  maximum 
debris  expansion,  and  at  =  80  [bottom  panel]  after  the  cloud  has  collapsed.  At  the 
earlier  time  most  of  the  debris  energy  has  gone  into  excluding  the  magnetic  field.  The  ions 
have  stopped  at  this  time  and  remain  fairly  cold.  When  the  magnetic  field  recollapses  and 
gives  its  energy  back  to  the  debris,  the  debris  ions  have  become  thermalized,  with  some 
ions  energized  to  about  AEo-  The  energy  distribution  function  can  be  characterized  at  this 
time  as  f{E)  ~  exp(— JEJ/JEJc),  with  Ec  cs  0.5J5<,. 

1.3.B  Scaling  Studies 

One  parameter  of  interest  in  VHANE  bursts  is  the  ratio  of  the  direct  ion  gyroradius  p, 
to  R3.  For  actual  VHANE’s  pi/Ri  -C  1,  but  numerous  experimental  studies,  including 
the  laser  experiments  at  NRL  [Ripin  et  al,  1987]  and  LLNL  [Dimonte  and  Wiley,  1991], 
and  various  simulations  [e.g.,  Winske,  1990]  have  considered  a  range  of  pi/R^  to  values 
exceeding  unity.  We  thus  consider  a  series  of  hybrid  simulations  in  the  manner  described 
above,  varying  the  value  of  Pi/Rz-  The  other  important  parameter  is  Md  =  yohAi  which 
we  will  consider  later. 

Figure  1.10  shows  the  debris  ions  in  the  x-y  plane  at  about  the  time  of  maximum  bubble 
size  for  four  veilues  of  PifR2-  At  small  Pi/R2  the  ions  have  compressed  into  a  thin  shell, 
and  flute  modes  have  begun  to  appear  on  the  siarface.  For  Pi/R2  =  0.63  the  instability, 
which  occurs  at  a  shorter  wavelength,  hais  already  heated  the  ions  in  the  shell  to  some 
extent  and  scattered  some  ions  to  larger  radius.  At  even  larger  pi/R2  =  1-25,  some  of 
the  ions  have  compressed  into  a  shell  that  shows  no  sign  of  instability,  while  more  ions 
have  been  scattered  to  an  even  larger  radius.  The  collapsed  debris  ions  at  later  time  are 
shown  in  Figure  1.11.  In  the  small  Pi/R2  cases  the  debris  is  confined  to  a  narrow  cylinder; 
at  intermediate  pi/R2  some  ions  remain  outside  the  main  cloud  in  a  thin  halo.  At  larger 
PilR2  a  significant  number  of  ions  do  not  recollapse  with  the  magnetic  field.  In  the  latter 
case,  however,  one  questions  the  validity  of  the  results.  In  the  full  particle  calciilations 
with  PilR2  >  1  a  strong  flute  instability  occurs  at  short  wavelengths  [Winske,  1990].  In 
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2.5i? 


— 2.5i?2  X  2.5/?2  ~  2.5i?2  X  2.5i?2 

Figure  1.6.  2-D  contour  plot  of  debris  density  at  various  times  for  an  expansion  at  Md  = 
0.5  with  Pi/R2  =  0.16,  showing  the  expansion  of  the  debris,  development  of  flute 
instability,  and  recollapse. 
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Figure  1.7.  3-D  perspective  plots  of  the  debris  density  and  debris  and  air  ion  energy  at 
Qit  =  80  for  the  same  calculation,  showing  highly  structured  nature  of  these  quantities. 
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Figure  1.8.  Radial  slice  of  the  debris  density  and  energy,  showing  debris  cloud  radius  ri 
and  halo  radius  r2. 
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Figure  1.10.  Debris  ions  in  x-y  space  at  about  the  time  of  maximum  bubble  size  for  four 
values  of  Pi/R2,  showing  some  ions  scattered  beyond  R2  when  Pi/R2  is  large. 
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Figxire  1.11.  Debris  ions  in  x-y  space  at  about  the  time  of  bubble  collapse  for  four  values 
of  PilR2-  At  large  pi/Ri  a  significant  nxunber  of  the  ions  do  not  recollapse  with  the 
main  debris  cloud. 
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this  regime  charge  separation  effects  are  significant.  The  quasi-neutral  hybrid  calculations 
presented  here,  however,  ignore  this  physics,  which  probably  explains  the  absence  of  an 
instability  in  this  case.  Nevertheless,  the  overall  tendency  for  some  ions  to  remain  at 
distances  on  the  order  of  1-2  pi  beyond  the  final  debris  column  size  is  evident. 

The  results  of  these  runs  are  summarized  in  Table  1.1.  In  order  to  vary  Pi/i?2  keeping 
R2  constant,  we  have  changed  the  mass  mo  and  the  initial  density  no  of  the  debris 
ions  compared  to  the  (fixed)  background  mo,  Mq.  Tabulated  for  these  four  runs  are  ri/i?2, 
T^Di/no,  EdxJEo,  and  EaiIEq  (radius  ri,  average  debris  density,  average  debris  ion  energy, 
and  average  ziir  ion  energy  for  r  <  rj),  as  well  as  the  corresponding  values,  r2,  nD2/”o, 
ED2/E0,  and  Ej\2/Eo,  for  ri  <  r  <  r2,  where  rig  is  the  background  density  and  Eg  is  the 
initial  debris  kinetic  energy  Eg  =  O.SmoVQ.  The  last  entry  in  the  table  is  /(r  >  ri),  the 
fraction  of  debris  ions  that  remain  outside  of  radius  rj.  Generally,  for  the  debris  in  the 
main  cloud  (r  <  rj),  ri/R2  decreases  with  Pi/R2,  Edi  is  approximately  constant,  and  npi 
and  Ej^i  increase  with  pi/R2-  In  the  hcdo  region  (ri  <  r  <  r2),  r2/R2,  n£)2/no,  EdiIEq, 
Ea2/Eo,  and  /(r  >  rj)  all  tend  to  increase  with  Pi/R2  (0.04  <  Pi/R2  <  0.63). 

Ignoring  the  somewhat  anomalous  case  with  Pi/R2  =  1-25,  one  would  conclude  that  almost 
all  of  the  debris  ions  tend  to  recollapse  with  the  magnetic  field  to  a  radius  r/  ~  0.5i22.  Less 
than  about  2%  of  the  debris  ions  remain  outside  the  main  cloud  in  the  VHANE  regime 
{pifR2  <  0.04).  Due  to  the  interaction  of  the  debris  with  the  magnetic  field  the  directed 
debris  energy  becomes  thermalized,  with  a  final  debris  ion  temperature  about  one-half  of 
the  initial  debris  kinetic  energy.  Air  ions  are  trapped  in  the  cloud  as  the  debris  collapses 
and  are  also  heated  to  a  comparable  temperature.  The  heating  of  the  air  ions  here  increases 
with  PilR2,  because  the  energy  is  shared  with  fewer  air  ions  relative  to  debris  ions.  The 
few  debris  (and  the  air)  ions  that  exist  in  the  halo  arovmd  the  debris  cloud  (halo  radius 
~  1  —  2pi)  2u:e  heated  to  similar  temperatirres.  In  an  actual  3-D  expansion  the  cloud  will 
zilso  elongate  along  the  magnetic  field,  leading  to  about  a  20%  small  final  radius  across  the 
magnetic  field  (i.e.,  r/  ~  0.4i?2)  when  the  field  collapses,  and  with  reduced  temperatures 
(the  energy  being  shared  by  all  three  degrees  of  freedom). 

1.4  Conclusions 

We  have  used  two  types  of  hybrid  codes  to  study  the  debris  in  a  VHANE  burst.  Czdcu- 
lations  with  a  new  hybrid  algorithm  confirm  our  previous  simiilations  results  and  most  of 
the  LLNL  and  BRA  calculations  that  in  a  3-D  expansion  debris  will  extend  out  to  about 
0.8  R3,  consistent  with  Kilb’s  estimate.  The  absence  or  presence  of  a  backgroimd  plasma 
makes  little  difference  in  the  expansion.  Two  dimensionzd  debris  expansions  perpendicular 
to  the  ambient  magnetic  field  have  been  carried  out  to  times  comparable  to  the  collapse  of 
the  magnetic  bubble.  The  debris  becomes  structured  in  the  field,  but  nevertheless  collapses 
to  a  final  radius,  ~  0.5i?2,  leaving  only  some  debris  ions  one  or  two  ion  gyroradii  beyond 
the  collapsed  cloud;  the  fraction  of  such  ions  decreases  as  pi/R2  is  reduced.  The  debris 
(and  air  ions  inside  the  debris  cloud)  become  thermalized,  with  energies  roughly  half  of  the 
directed  debris  kinetic  energy.  The  ratio  of  pijRz  has  little  overall  effect  on  the  thermal- 
ization  or  the  development  of  the  flute  instability,  as  long  sis  PilR2  <  1-  This  result  seems 
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Table  1.1  Debris  zind  Air  Characteristics  for  Hybrid  Runs  at  Various  Pi/R2- 


Run 

1 

2 

3 

4 

nofrio 

163.0 

41.0 

10.2 

5.1 

^dI^o 

1.0 

4.0 

16.0 

32.0 

PijRa 

0.04 

0.16 

0.63 

1.3 

ri/RB 

0.61 

0.53 

0.45 

0.25 

noi/T^o 

0.08 

0.09 

0.14 

0.22 

Edi/Eo 

0.63 

0.51 

0.55 

0.62 

EaI  / Eg 

1.39 

2.30 

4.32 

2.71 

r2lRB 

0.19 

0.22 

0.45 

2.34 

^02(^0 

0.002 

0.003 

0.009 

0.0003 

Ed2  /  Eg 

0.34 

0.69 

0.91 

0.27 

Ea2/Eo 

1.39 

2.30 

4.32 

2.71 

/(r  >  n) 

0.017 

0.032 

0.16 

0.48 
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somewhat  contrary  to  our  earlier  results  obtained  with  full  particle  codes  [Winske,  1990], 
However,  in  the  particle  calculations  short  wavelength  (lower  hybrid)  modes  dominated, 
which  are  suppressed  as  pi/R2  decreases.  The  hybrid  simulations,  on  the  other  hand, 
involve  longer  wavelength  modes  that  appau"ently  are  less  susceptible  to  Pi/R2  effects. 

While  the  hybrid  results  presented  here  have  been  done  only  in  2-D,  one  expects  that  in 
3-D  the  final  debris  radius  to  be  slighly  smaller,  due  to  the  elongation  of  the  cloud  zdong 
B,  yielding  ~  0.4  *  O.SR3  ~  0.3/?3,  which  interestingly  is  similar  to  that  which  Brecht 
[private  communication]  gets  in  some  cases  when  the  debris  anomjdously  slips  through  the 
field  without  making  much  of  a  magnetic  cavity.  And  in  3-D  one  expects  the  debris  to 
thermalize  in  all  three  directions  to  a  final  energy  of  about  l/3Eo  (rather  than  l/2Eo  in 
2-D).  Yet  to  be  determined  is  the  slope  of  the  more  energetic  part  of  the  debris  energy 
spectrum  and  the  eventual  deposition  of  the  debris  energy  in  the  ionosphere. 
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2.  Energy  Deposition  Spectrum  in  a  HANE  Burst 
2.1  Introduction 

The  essential  question  for  nuclear  explosions  at  high  adtitudes  is  how  the  weapon  energy 
is  deposited  in  the  atmosphere  when  the  ambient  air  density  is  so  low  that,  when  fully 
ionized,  the  medium  is  collisionless.  Much  of  the  energy  goes  off  as  ionizing  radiation,  but 
some  fraction  remains  as  kinetic  energy  in  the  debris.  The  initial  expansion  velocity  of 
the  debris  ions  is  much  greater  than  the  local  Alfven  speed  so  that  a  high  Mach  number 
collisionless  shock  forms,  which  causes  the  debris  to  stop  as  it  gives  up  its  energy  to  form 
the  shock,  and  energizes  the  air  ions  that  are  eventually  deposited  in  the  atmosphere.  This 
complex  process  is  modeled  differently  in  various  “first  principles”  codes.  In  CMHD  [Kilb 
and  rJlen,  1978]  niimerical  viscosity  is  included  to  stably  model  the  air-shock  process  and 
through  this  viscous  interaction  the  air  ions  are  thermalized,  resulting  in  a  rather  soft 
energy  spectrum.  In  the  BRA  SHYPS  code  [Brecht  and  Thomas,  1988]  the  ion  kinetics  of 
the  debris-air-shock  are  modeled  self-consistently,  resulting  in  a  nonclassical  collisionless 
shock  and  a  more  energetic  deposition  spectrum. 

Recently,  the  consequences  of  these  two  different  energy  deposition  spectra  have  been 
compared  using  the  systems  code  NORSE  [Thompson,  1991].  The  spectra  used  as  initial 
conditions  in  these  two  calculations  differ  in  a  number  of  aspects:  total  energy  deposited 
(CMHD  deposits  about  three  times  as  much  energy  as  SHYPS),  radial  distribution  of  the 
energy  (SHYPS  deposits  more  energy  further  from  the  burst  point),  and  velocity  distri¬ 
bution  (while  the  CMHD  spectrum  falls  off  monotonically,  the  SHYPS  spectrum  has  two 
peaks  and  extends  to  higher  velocity).  Not  suprisingly,  the  two  deposition  models  give 
rise  to  significantly  different  nuclear  phenomenologies:  different  ionization  levels,  electron 
temperatures,  and  atmospheric  heave.  Of  particular  interest  here  are  the  resulting  differ¬ 
ences  in  the  electron  density,  which  at  certain  altitudes  and  times  are  about  two  orders 
of  magnitude  larger  in  the  calculations  using  the  CMHD  spectra.  In  part  this  is  due  to 
the  greater  total  energy  deposited  in  the  CMHD  run  and  its  concentration  closer  to  the 
burst  point  field  line.  It  is  also  due  to  the  softer  CMHD  spectrum  that  deposits  the  energy 
higher  in  the  atmosphere  and  shares  it  with  fewer  background  ions,  yielding  higher  electron 
temperatures. 

In  order  to  help  understand  these  differences  and  perhaps  reconcile  them,  we  were  eisked 
by  DNA  to  provide  an  independent  verification  of  the  energy  deposition  spectra  in  HANE 
events.  Rather  than  reproduce  the  BRA  result,  we  have  carried  out  more  restricted  kinetic 
calculations  (which  we  coiild  do  on  a  rather  short  time  scale)  to  demonstrate  the  essential 
features  of  the  spectrum.  These  simulations  are  described  in  Section  2.2.  These  features 
are  also  corroborated  by  observations  of  high  Mach  number  collisionless  shocks  in  space; 
we  briefly  review  some  of  the  major  findings  in  Section  2.3.  Finally,  in  Section  2.4  we 
summarize  our  results.  We  also  describe  other  effects  that  may  tend  to  bring  the  two 
models  into  better  agreement  and  discuss  some  further  tests  that  can  be  done  in  the 
systems  codes  to  understzind  the  differences. 
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2.2  Simulations 


Our  purpose  here  is  to  show  the  origin  of  the  energized  air  spectrum,  which  is  not  a 
deposition  spectrum  per  ae,  but  which  contains  the  essentieJ  features.  As  such  we  carry 
out  the  calculations  in  restricted  geometries.  The  simulations  are  done  with  a  hybrid  code 
that  treats  the  ions  kinetically  by  mezms  of  paxticle-in-cell  techniques  zmd  the  electrons 
35  a  massless  fluid,  similar  to  that  used  by  Clark  et  al.  [1973],  Goodrich  ei  al.  [1985], 
and  Thomas  and  Brecht  [1986].  We  describe  simple  planar  and  cyiindriczd  calculations  in 
one  and  two  spatial  dimensions.  In  the  planzir  calculations  [Winske  and  Quest,  1986],  the 
debris  initially  is  uniform  over  the  region  0  <  x  <  Ld  and  flows  in  the  +x  direction  with 
initial  velocity  Vp.  The  air  occupies  the  remaining  region  Lq  <  x  <  L,  is  uniform,  and 
at  rest  initially.  A  constant  magnetic  field  of  magnitude  Bo,  which  lies  m  the  x-y  plane, 
fills  the  region.  A  closed  system  with  reflecting  particle  boimdaries,  on  which  the  electric 
field  is  zero,  is  assumed.  In  the  two-dimensional  calculations  similar  initial  and  boundary 
conditions  along  x  are  employed;  in  the  y  direction  the  plasma  is  uniform  and  periodic 
boundary  conditions  axe  imposed. 

We  begin  with  a  simple  1-D  simulation  that  shows  the  essential  features  of  the  interaction. 
In  this  case  the  debris  Mach  number,  with  the  Alfven  speed  defined  with  the  air  parameters 
iv\  =  BljATTriAmA),  is  16.  The  system  length  is  L  =  —  47rnAe^/m^],  the 

cell  size  8  —  0.25c/u;yi,  rn£)lmA  =  2,  £d  =  20c/u)a,  I^d  =  =  /^  =  T  c/va  = 

^aI^a  =  4000,  and  the  resistivity  rj  =  10~®u;^/47r.  Figure  2.1  shows  the  normalized 
debris  density  no(x)/nDit  =  0)  (dashed  line),  air  density  n4{x)/nx(f  =  0)  (sohd  line), 
and  magnetic  field  magnitude  B{x)/Bo  (dotted  line)  as  a  fimction  of  x.  (Note  that  the 
normalized  B  field  is  given  on  a  different  scale  so  that  it  can  be  more  easily  seen  in  the 
figures.)  Frcunes  at  three  different  times  are  displayed:  =  0,flAt  =  2.5,  and  =  10 

{Q.A  =  eBo/m.4C=gyrofrequency  of  air  ions).  The  =  0  show  the  initied  configuration. 
The  steep  gradients  give  rise  to  a  small  pulse  on  the  magnetic  field,  which  is  not  easily 
seen  in  the  figure.  Other  initial  configurations  with  smoother  density  profiles  or  compressed 
magnetic  fields  in  the  debris  region  have  also  been  used.  Except  for  some  slight  differences 
at  very  early  time,  which  zue  to  be  expected,  the  long  time  behavior  is  similar  [e.g.,  Clark  et 
al,  1973].  By  =  2.5  the  debris,  the  air,  and  the  magnetic  field  have  all  been  strongly 
compressed,  with  a  shock-like  structure  already  formed.  At  later  times,  =  10,  the 
shock  propagates  to  the  right  through  the  unperturbed  air,  while  the  debris  spreads  out 
and  eventually  stops.  Note  that  in  this  strictly  perpendicular  configuration  that  there  is 
little  mixing  of  the  debris  and  air  plasmas;  a  quantitative  mezisure  of  which  is  given  in 
Wmske  and  Quest  [1986].  Both  the  magnetic  field  zmd  air  density  show  large  pulses,  czdled 
“overshoots”,  at  the  leading  edge  of  the  shock. 

The  dynzimics  of  the  interaction  is  better  seen  in  —  x  phase  space  for  both  ion  species 
at  these  same  times,  as  given  in  Figure  2.2  Again,  at  <  =  0  both  species  have  uniform 
densities  and  velocities.  At  Qa^  =  2.5,  the  compression  of  the  magnetic  field  causes  the  air 
ions  at  the  debris-air  interface  to  be  picked  up,  gyrate  in  the  compressed  field  to  velocities 
larger  than  Vo,  while  the  debris  motion  is  somewhat  more  complicated  at  this  early  time. 
The  dynamics  at  very  early  times  have  been  described  by  Goodrich  et  al  [1985].  At  later 
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Figure  2.1.  Results  of  1-D  planar  simulations  at  varioiis  times  showing  air  density  (solid 
curves),  debris  density  (dashed  curves),  and  magnetic  field  magnitude  (dotted  curves) 
normalized  to  their  initial  values  for  Mn  =  16,  <5  =  1,  ^  =  90*. 


times,  a  fraction  of  the  backgroimd  ions  interacting  with  the  leading  pulse  are  strongly 
scattered.  This  ion  reflection  process  generally  occurs  at  supercritical,  quasiperpendicular 
collisionless  shocks  and  has  been  well  described  by  Leroy  et  al.  [1981],  Goodrich  [1985],  and 
for  high  Mach  number  shocks  by  Quest  [1986].  These  reflected  ions  determine  the  main 
features  of  the  shock  structure  and  give  rise  to  the  harder  SHYPS  deposition  spectrum. 
The  debris,  on  the  other  hand,  has  essentially  stopped.  Some  debris  ions  still  have  positive 
X  velocities,  but  other  have  been  reflected  at  the  front  edge  of  the  debris  and  have  been 
heated. 

The  process  of  shock  formation  and  energization  of  the  air  ions  is,  of  course,  not  just 
restricted  to  planar  calculations.  Figures  2.3  and  2.4  show  the  results  of  a  cylindrical 
expansion  of  debris  ions  into  background  eiir,  again  with  an  initial  Mach  number  of  16. 
Figure  2.3  shows  the  magnetic  field  (directed  out  of  the  plane  of  the  simulation)  about 
the  time  the  debris  (initially  a  dense  cyhnder  expanding  radially)  has  come  to  rest,  with  a 
cylindrical  shock  expanding  outward.  In  Figure  2.4  radial  cuts  of  B  and  the  background 
ions  in  Vi  —  a;  and  Vy  —  x  phase  space  at  the  same  time  [left  panels]  show  the  shock  fronts  at 
at  I  ~  10  and  x  ~  90,  with  corresponding  energetized  background  ions.  The  right  panels 
of  Figure  2.4  show  the  air  and  debris  ions  in  the  x-y  plane  at  the  same  time.  The  debris 
ions  have  come  to  rest,  while  most  of  the  background  ions  have  been  swept  out  by  the 
shock.  Outside  the  shock,  the  backgroimd  ions  are  still  undisturbed. 

Planar  simulations  with  two  spatial  dimensions  have  also  been  run  using  the  seune  initial 
conditions.  These  calculations  show  the  same  ion  reflection  and  energization  process  as  the 
1-D  simulations  and  because  they  include  the  magnetic  field  in  the  plane  of  the  singulation 
{B  —  BoV),  they  further  show  how  ions  axe  scattered  in  the  direction  of  the  magnetic  field 
by  an  electromagnetic  instability.  These  energetized  air  ions  can  then  stream  along  the 
magnetic  field  and  be  deposited  at  lower  altitudes.  As  with  the  first  example  presented  in 
Figures  2.1  and  2.2,  the  2-D  example  is  characterized  by  Mp  =  16,/?  =  1,^  =  90°  with 
Ld  =  ^c/uA  and  no  In  a  =  2.  The  system  size  is  200  x  32  with  the  cell  size  Ai  =  Aj,  = 
Q.^c/ua-  Results  at  =  7.5  axe  displayed  in  the  figmes.  Figure  2.5  shows  —  x  and 
Vy  —  x  phaise  space  for  both  debris  ions  [left  panels]  and  edr  ions  [right  panels].  In  the  air  ion 
phase  space  one  again  sees  a  well  formed  shock;  the  reflection  process  is  somewhat  less  clear 
because  the  figure  averages  over  the  entire  y  direction.  Most  significantly  is  that  heating  of 
the  air  ions  in  Vy  along  the  magnetic  field  is  also  seen.  This  effect  is  due  to  the  Alfven  ion 
cyclotron  instability,  as  such  an  effect  does  not  occur  in  the  one-dimensioneJ  simulations 
[Tanaka  et  al,  1983;  Thomas  and  Brecht,  1986].  The  heating  occurs  behind  the  shock  front, 
and  leads  to  parallel  velocities  not  much  different  that  those  in  the  perpendicular  direction. 
The  debris,  on  the  other  hand,  shows  little  parallel  heating,  which  is  not  suprising,  as  the 
Alfven  ion  cyclotron  instability  occurs  only  near  the  front  of  the  shock.  Profiles  in  x 
(averaged  over  all  y)  of  the  normzilized  debris  (detshed  line)  and  air  (solid  line)  [top  panel] 
and  magnetic  field  magnitude  [bottom  panel]  are  given  in  Figure  2.6  at  the  same  time.  In 
this  case  there  is  slightly  more  debris-air  mixing  with  some  eur  ions  left  behind.  The  air 
density  and  magnetic  field  show  a  well  formed  shock,  similar  to  Figure  2.1. 
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Figure  2.4,  [left  panels]  Cuts  through  center  of  simulation  region  showing  magnetic  field 
profile  and  background  ion  v*  —  z  and  —  x  phase  space;  [right  panels]  debris  and 
backgroimd  ions  in  x-y  at  the  same  time  of  the  calculation  shown  in  Figure  2.3. 
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Figure  2.5.  Results  of  a  planar  2-D  simxilation  with  parameters:  Md  =  16,0  =  1,  Ld  = 
10c/u)A,  B  =  90"  at  flyit  =  7.5,  showing  Vg  —  x  and  Vp  —  x  phase  space  for  debris  ions 
(left  panels)  and  air  ions  (right  panels). 
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Contour  plots  of  the  debris  and  zur  densities  are  displayed  in  Figure  2.7.  One  notes  im¬ 
mediately  that  the  shock  front  is  slightly  rippled.  In  addition,  the  air  density  shows  some 
well  defined  wave  structure  just  behind  the  shock,  with  a  peak  density  of  11.  The  con¬ 
tours,  which  are  not  labelled,  rim  from  0.0  to  1.6  in  steps  of  0.1  for  the  debris  and  0.0 
to  11.0  in  steps  of  1.0  for  the  air.  They  indicate  more  clearly  the  lack  of  mixing  of  the 
two  plasmas  and  the  absence  of  any  structure  in  the  debris  ions.  Corresponding  2-D  plots 
of  the  magnetic  field  components  By  and  Bz  are  given  in  Figure  2.8.  Contours  for  By 
range  from  0.0  to  5.0  in  steps  of  0.5,  while  those  for  range  from  -2.0  to  2.0  with  similar 
increments.  While  the  turbulence  in  By  is  less  evident  because  of  the  macroscopic  shock 
profile,  the  level  of  fluctuations  in  By  is  comparable  to  that  in  Bz  (as  well  as  Bx),  with 
SB /Bo  ~  2.0.  The  primary  area  of  wave  activity,  as  best  seen  in  Bz,  is  slightly  behind 
the  shock  front,  the  region  where  the  reflected  ions  become  isotropized  by  the  Alfven  ion 
cyclotron  instability. 

The  important  new  results  of  the  calculation  aure  shown  in  Figure  2.9.  From  the  2-D 
simulation  shown  in  Figures  2. 5-2.8  we  calculate  the  parallel  velocity  spectrum  f{v\\)  [top 
panel]  and  the  related  “deposition  spectrum”  dE/dv\\  [bottom  panel]  of  the  air  ions  just 
behind  the  shock  front,  as  a  function  of  U||  at  the  same  time  Q,it  =  7.5.  The  spectra 
are  normadized  to  the  shock  velocity  u,  8va)-  The  solid  curve  in  the  top  panel  is 
the  corresponding  Maocwellian  distribution.  At  low  U||,  /(u||)  is  Maxwellian,  but  deviates 
significantly  at  high  U||  due  to  the  energized  reflected  ions.  Because  of  the  high  energy  tail, 
dE/dv\\  is  double  peaked  (compare  with  model  CMHD,  NHL,  and  BRA  spectra  shown 
in  Figure  2.10,  from  Brecht  et  al.  [1990]).  As  is  readily  seen,  the  spectrmn  here  agree 
quantitatively  with  the  SHYPS  spectra  in  two  important  aspects.  First,  the  spectrum 
generally  has  two  peaks:  one  at  low  energy,  due  to  the  thermal  ions,  and  one  at  higher 
energy,  due  to  the  reflected  ions.  Secondly,  the  spectrum  extends  out  to  about  2t;,  (which 
in  this  caise  is  about  the  initial  debris  velocity).  (Had  we  driven  the  shock  harder,  the  shock 
velocity  would  have  been  about  the  debris  velocity  end  the  spectrum  would  then  extend 
to  about  twice  the  initial  debris  velocity,  as  it  does  in  the  SHYPS  calculation.)  These  two 
features  are  generic  properties  of  the  kinetic  apects  of  the  interaction  of  the  air  ions  with 
the  shock  [e.g.,  Tanaka  tt  al,  1983]  and  are  not  affected  much  by  the  limited  geometry  of 
the  present  calculation. 

That  the  spectrum  we  find  agrees  with  the  BRA  calculations  is  not  imexpected,  as  the 
two  simulation  codes  contain  the  same  physics,  even  though  we  have  done  a  planar  shock 
2ind  BRA  has  carried  out  a  more  complete  calculation.  Also,  one  should  note  that  our 
“deposition  spectrum”  is  computed  is  the  source  region,  not  the  deposition  region.  One 
expects  fewer  mr  ions  at  lower  un  to  reach  the  deposition  region  early  in  the  explosion 
and  hence  there  will  be  an  absence  of  ions  at  small  dE/dv\\,  enhancing  the  double  peaked 
structure  of  the  spectrum.  Also,  with  our  planar  shock,  the  shock  does  not  slow  down  very 
much  in  time,  so  that  our  spectrum  does  not  evolve  in  time  as  does  the  BRA  spectrum.  In 
a  more  realistic  calculation,  as  the  shock  expands,  it  slows  down.  At  lower  Mach  number, 
the  number  of  air  ions  reflected  by  the  shock  front  decreases,  resulting  in  a  weaJcer  peak 
at  higher  energy  in  the  spectrum. 
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Figure  2.7.  Contour  plots  in  x  and  y  of  air  [top]  and  debris  [bottom]  ion  densities  normalized 
to  their  initial  values;  contoxus  arc  evenly  spaces  and  vary  as  0.0,0.05,0.1, ,  ,0.6  for 
debris  and  0.0, 0.5, 1.0, , ,  7.0  for  air. 
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Figtire  2.8.  Contour  plots  in  x  and  y  of  B,  [top]  and  B,  [bottom]  normalized  to  Bo; 
contoiirs  arc  evenly  spaces  and  vary  as  0.0, 0.5, 1.0, ,  ,5.0  for  B^  and  —2.0,  —1.5, ,  ,2.0 
for  B,. 
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Figure  2.9.  [top  panel]  Parallel  velocity  spectnnn,  /(ny),  for  the  air  ions  just  behind 
the  shock;  [bottom  panel]  deposition  spectrum,  dEjdv^y  at  the  same  time  for  the 
simulation  in  Figures  2.5-2.8.  Both  spectra  are  plotted  versus  uy,  normalized  to  the 
shock  speed,  v,. 
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Figure  2.10.  Model  deposition  spectra  from  CMHD  and  NRL  models  and  computed  spectra 
from  SHYPS  (from  Brecht  ei  oi.,  1990). 
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2.3  Evidence  from  Space  Data 

The  above  picture  of  a  rather  haxd  deposition  spectrum  resulting  from  the  interaction  of 
an  outgoing  magnetic  shock  with  background  air  ions  is  supported  by  data  from  various 
spacecraft.  We  briefly  show  three  examples  that  corroborate  the  results  presented  previ¬ 
ously.  Detailed  description  of  the  results  and  supporting  modeling  efforts  are  found  in  the 
references. 

1.  The  collisionless  shock  that  results  from  a  high  Mach  number  debris  expansion  is  non- 
laminar.  Unlike  a  classical  MHD  shock,  the  magnetic  field  profile  is  chziracterized  by  a 
large  “overshoot”  at  the  leading  edge  of  the  shock.  This  overshoot  is  seen  in  numerous 
crossings  of  the  bow  shock  of  the  Earth  and  other  planets,  especially  the  outer  planets 
where  the  Mach  number  is  higher.  Figure  2.11  is  from  Bagenal  ei  al.  [1987],  showing  the 
bow  shock  of  Uranus.  The  Mach  number  here  is  23.  The  shock  front  exhibits  a  very  large 
overshoot;  behind  the  shock  front,  the  magnetic  field  reduces  to  the  level  expected  from 
the  Rankine-Hugoniot  relations.  These  results  are  also  relevant  to  a  recent  einalysis  of  the 
STARFISH  magnetic  field  measurements  [Brecht  and  Thomas,  1988). 

2.  The  interaction  of  the  shock  produces  two  populations  of  air  ions.  One  population  is 
transmitted  through  the  shock  and  heated.  The  other  population  is  reflected  and  ener¬ 
gized.  These  two  populations  then  give  rise  to  the  double  peaked  nature  of  the  deposition 
spectrum.  Spacecraft  observations  [Sckopke  et  o/.,  1983]  show  this  reflection  process  and 
the  two  distinct  populations  in  great  detail  [Figure  2.12].  Simulations  of  this  process 
demonstrate  the  dynamics  quite  clearly  [Leroy  et  al,  1981]. 

3.  The  subsequent  thermaJization  of  the  reflected  ions  and  their  energization  along  the 
magnetic  field  by  the  Alfven  ion  cyclotron  instability  produce  a  highly  nonthermal  distri¬ 
bution  which  persists  far  downstream  of  the  shock  (Figure  2.13).  Again,  the  observations 
of  Sckopke  ei  al  [1983]  show  the  non-Maxwellian  natiure  of  the  ion  distribution  and  the 
persistence  of  an  energized  ion  component  in  the  downstream  region,  which  is  the  origin 
of  the  harder  SHYPS  energy  spectrum.  Simulations  confirm  both  the  spectrum  shape  and 
the  persistence  of  this  effect  in  the  downstream  [Winske  and  Quest,  1988]. 

2.4  Conclusions 

We  have  used  hybrid  simulations  of  high  Mach  number  debris-air  ion  interactions  to  verify 
the  essental  features  of  the  SHYPS  energy  deposition  spectrum  that  have  been  modeled 
in  NORSE  calculations  by  Thompson  [1991].  Contrary  to  the  CMHD  model,  the  energy 
deposition  spectra  are  characterized  by  two  peaks  at  finite  V||  and  a  maximum  value  of  vn 
that  extends  to  above  the  initial  debris  velocity.  We  have  also  shown  examples  of  space 
data  that  lend  additional,  independent  support  to  this  picture. 

While  the  NORSE  calculations  are  a  good  start  in  putting  some  of  the  kinetic  aspects 
of  debris-air  interactions  into  systems  codes  and  showing  that  significant  differences  are 
possible,  we  believe  that  these  effects  should  be  pursued  in  parallel  with  other  modeling 
efforts  enumerated  below. 


36 


0700  0800  0900 

SCET 


Fig.  S.  The  9.2-1  imigei  of  the  mignctic  Seld  migniadf  (8),  directiao  (aglet  X,  8),  nd  nnt  in  heliognphic 
co-ordinttes.  The  angle  X  ii  defined  by  taX  ■  BjtBt.  The  ogle  8  is  defined  by  tin8  ■  BtJB.  Tlie  vaiuei  cX  the 
parameten  predicted  by  the  Rankine-Hugaaiot  (R-H)  reUtiont  are  shown  by  hoiizonul  lines. 


Figure  2.11.  Profiles  of  the  magnetic  field  magnitude  and  other  related  quantities  at  the 
crossing  of  the  Uranus  bow  shock,  showing  large  overshoot  [from  Bagenal  ti  al.,  1987]. 
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Fig.  2.  ion  velocity  iptoe  diitributioBi  for  tbc  inbound  thoefc  croaing  on  November  7, 1977,  spanning  the  range  from 
the  pure  solar  wind  to  the  shock  ramp.  To  emphasize  the  evolution  of  the  distributions,  only  every  other  of  those 
measured  is  shovm.  The  relative  positions  where  the  measurements  were  made  arc  indicated  by  the  dots  on  the  plasma 
density  profiles  shown  as  insets  ^at  cover  an  interval  of  2.5  min.  The  distributions  are  shown  u  contours  of  constant 
phase  space  density  in  two^imensional  velocity  space,  /(e.  p)  •  const,  where  c  »  the  ion  speed  and  p  is  the  ecliptic 
azimuth  of  the  velocity  vector,  with  p  -  0  pointing  toward  the  sun  on  the  left.  Density  levels  represented  by  adjacent 
contours  differ  by  an  order  of  magnitude.  The  plus  symbols  indicstc  the  coordinate  origin  c  >  0  in  spacecraft  (or, 
approsimately,  the  shock)  franw  of  reference  with  the  vdocity  scale  given  in  the  first  (fame.  The  dashed  line  in  this  frame  is 
the  projection  of  a  vector  parallel  to  the  shock  front  and  perpendicular  to  the  magnetic  field,  onto  the  ecliptic  plane.  The 
asterisks  in  frames  2  and  9  indicate  the  spaed  and  azimuth  of  specularly  refracted  ions  at  their  forward  turning  point  and 
near  the  sbo^  ramp,  respectively. 


Figtire  2.12.  Ion  phase  space  measured  at  the  Earth’s  bow  shock  showing  two  ion  popula¬ 
tions  [from  Sekopke  ti  ai,  1983]. 
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Fig.  IS.  Time  sequence  of  one-dimensional  ion  distributions 
downstream  of  the  November  7, 1977,  shodc,  beginning  a  few  seconds 
after  the  shock  crossing.  The  distributions  are  shown  as  cuts,  /(£, 
^  «  const),  along  the  direction  O  of  the  bulk  flow,  with  Eccv*  given 
in  keV,  and  /  in  units  of  cm*  s'*.  Secondary  mMima  (rf  /  due  to 
gyrating  ions  are  present  only  for  a  short  time,  but  a  broad  shoulder 
of  /  persists,  nearly  undianged,  for  more  than  one  hour.  The  dashed 
curve  represents  a  Maxwellian  fit  to  the  low-energy  part  of  the  distri¬ 
bution. 


Figure  2.13.  Ion  energy  distributions  downstream  of  a  bow  shock  crossing  showing  persis¬ 
tence  of  non-Maxwellian,  high  energy  ion  component  [from  Sckopke  et  al,  1983]. 
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1.  The  use  of  direct  output  at  discrete  times  from  SHYPS  to  initialize  NORSE  has  been 
an  expedient  way  to  use  the  SHYPS  data.  A  more  useful  format  would  be  einalytic  expres¬ 
sions  for  the  energy  deposition  spectra  that  are  derived  from  the  properties  of  the  air-shock 
interaction.  This  would  then  allow  a  systematic  study  of  kinetic  effects  on  nuclear  phe¬ 
nomenology  based  on  the  various  system  codes  for  a  range  of  burst  parameters.  This 
so-called  “turning  model”  could  also  be  used  in  CMHD  to  distribute  the  energy  resulting 
from  viscous  heating  to  the  air  ions. 

2.  One  way  to  reconcile  the  differences  between  the  two  models  might  be  to  try  the 
deposition  model  of  Witt  and  Kilb  [1988]  that  uses  blast  wave  formulas  to  calulate  shock 
parameters  zuid  the  Rankine-Hugoniot  relations  to  distribute  the  energy.  These  formulas 
could  be  generalized  to  introduce  a  population  of  energetic  reflected  ions,  consistent  with 
the  hybrid  simulations. 

3.  One  effect  that  will  bring  the  two  models  into  better  agreement  is  the  possibility  of 
futher  electromagnetic  interactions  of  the  energized  air  ions  as  they  propagate  down  the 
field  lines.  At  the  debris  modeling  meeting  last  fall  we  showed  preliminary  calculations  of 
electromagnetic  ion/ion  instabilities  that  could  slow  and  thermalize  air  ions,  leading  to  a 
softer  energy  spectrum  and  deposition  at  higher  altitudes  due  to  such  anomalous  effects. 
Some  of  this  work  will  be  presented  later  in  this  report. 

4.  The  initial  NORSE  comparisons  using  the  SHYPS  and  CMHD  spectra  contain  a  number 
of  important  differences:  (a)  total  energy  deposited,  (b)  energy  spectrum  shape,  and  (c) 
energy  deposited  with  distance  from  the  burst  point.  All  three  of  these  differences  can 
affect  the  resulting  phenomena  in  various  ways.  It  would  be  useful  to  carry  out  the  systems 
calculations  varying  only  one  effect  at  a  time  to  see  how  these  processes  result  in  which 
differences  to  the  phenomenology  and  how  sensitive  those  differences  are  to  the  models. 

5.  As  noted  by  Thompson  [1991],  the  present  NORSE  calculations  do  not  include  debris 
energy  that  remains  imdeposited,  which  for  the  SHYPS  calculations  czm  be  significzint. 
What  is  needed  is  a  model  for  where  the  debris  is  and  what  its  energy  characteristics  are 
[as  heis  become  very  apparent  at  the  last  two  debris  modeling  meetings],  zis  well  as  the 
capability  to  include  such  undeposited  debris  energy  in  the  systems  codes. 

6.  Finally,  the  study  of  high  Mach  number  debris  air  interactions  in  the  collisonless  limit 
using  the  AGNES  facility  would  be  highly  desirable.  Properly  designed  experiments  could 
measure  the  deposition  spectra  and  provide  a  controllable  environment  for  doing  detailed 
modeling  studies. 
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3.  Debris/ Air  Interactions 

3.1  Introduction 

The  collisionless  interaction  of  energetic  debris  ions  with  ionized  etir  ions  in  HANE  bursts 
remziins  an  important  problem.  While  numerous  aspects  of  this  interaction  have  been 
studied  for  more  than  two  decades,  improvements  in  numerical  simulation  methods  and  of 
understanding  of  the  underlying  physical  processes  have  led  to  continuing  investigations. 
Indeed,  Thompson’s  [1991]  recent  comparisons  of  the  effect  of  different  eaurly  time  depo¬ 
sition  spectra  in  NORSE  calculations  on  late  time  characteristics  of  electron  temperature 
and  density  indicate  significant  differences  between  two  early  time  models  and  a  continued 
need  for  refinement  of  early  time  calculations.  In  this  section  we  concentrate  on  debris 
characteristics,  such  as  spatial  distribution,  energy  spectra,  ard  interactions  with  air  ions, 
as  obtained  from  hybrid  (particle  ion,  massless  fluid  electron)  simulations.  These  calcxila- 
tions  are  relatively  simple,  in  that  they  involve  only  two  spatial  dimensions  and  expansion 
perpendicular  to  a  vmiform  geomagnetic  field.  The  plan  here  is  to  look  for  scaling  relations 
as  a  function  of  initial  debris  Mach  number.  Also,  the  simulations  zire  used  to  develop 
new  diagnostics  to  characterize  the  properties  of  the  debris.  Later,  we  plan  to  csmry  out 
more  extensive  simulations  to  test  the  scaling  relations  and  include  additioned  physics.  W’e 
will  also  compare  these  results  with  simulation  studies  of  Brecht  [unpublished,  1991]  and 
Simonson  and  Hewett  [1991]. 

3.2  Simulations 

Figure  3.1  shows  the  basic  configuration  used  in  the  simulations.  A  debris  cloud  with  a 
Gaussian  density  profile  [nD(r)  =  noo  exp(— r^/a^),  a  ~  2c/u;i]  is  situated  in  the  center 
of  the  simulation  domain  (0  <  r  <  Lx]  0  <  y  <  Ly]  with  Lx  =  Ly  =  lOOc/wj).  The 
debris  ions  with  mass  mo  expand  radially  with  velocity  Vq  and  have  a  thermaJ  velocity 
VtD^  expressed  in  terms  of  an  ion  beta  The  debris  expands  into  a 

uniform  background  plasma,  with  density  and  ion  mass  rrio.  An  ambient  magnetic  field 
BgZ  is  perpendicular  to  the  simulation  plane.  With  these  backgroimd  parameters  we  define 
the  ion  inertial  length,  c/w^,  where  c  is  the  speed  of  light  and  is  the  ion  plasma  frequency 
(=  [47rnoC^/mo]^^^),  and  the  Alfven  speed  va  =  J3o(47rnom<,)“^^^.  Time  is  measured  in 
terms  of  the  inverse  of  the  ion  cyclotron  frequency,  =  eBoIrrioC.  In  these  calculations 
we  take  mo  =  2mo,  —  0.01,  while  the  ratio  of  udoI^o  is  varied  (see  Table  3.1). 

At  the  boundary  of  the  simulation  the  electric  field  is  set  equal  to  zero,  as  are  the  normal 
derivatives  of  the  magnetic  field  components.  This  allows  electromagnetic  (magnetosonic) 
waves  to  be  partially  reflected  at  the  boundaries.  An  alternative  plem  would  be  to  let  aJl 
the  waves  out  of  the  system,  but  the  boundary  conditions  that  we  use  have  the  advantage 
that  the  recovery  of  the  magnetic  field  is  modeled,  at  least  for  the  lower  MaxJi  number 
expzmsions.  As  will  be  seen,  in  the  higher  Mach  number  cases,  these  boimdary  conditions 
give  rise  to  a  reflected  shock  wave.  This  is  not  quite  appropriate  for  HANEs,  although  it 
is  an  interesting  way  to  model  multi-bursts.  To  mitigate  the  reflection  effect,  the  normal 
velocities  of  energized  background  zur  ions  are  set  equzd  to  zero  at  the  boimdary.  For  our 
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Table  3.1.  Parameters  and  results  of  simulations  at  different  debris  Mach  number. 


Md 

0.2 

0.4 

0.75 

1.0 

2.0 

4.0 

8.0 

^Dof^o 

512 

128 

56.9 

32 

32 

32 

32 

Rsicl^i) 

12.8 

12.8 

16 

16 

32 

64 

128 

RM{c/uJi) 

64 

32 

21.3 

16 

16 

16 

16 

Eo/0.bmov\ 

4 

16 

56.3 

100 

400 

1600 

6400 

r\/r. 

0.43 

0.49 

0.47 

0.47 

0.91 

1.06 

1.05 

0.37 

0.30 

0.27 

0.30 

0.20 

0.24 

0.92 

”Dl/”Do 

0.13 

0.12 

0.08 

0.07 

0.02 

0.003 

0.005 

T^D2/«Do 

0.006 

0.005 

0.01 

0.005 

0.001 

0.001 

0.003 

riAl/rio 

0.55 

0.66 

0.77 

0.60 

0.31 

0.10 

0.035 

^A2l^o 

0.98 

1.05 

0.96 

1.13 

0.91 

1.07 

0.90 

EdiIRo 

1.65 

0.75 

0.29 

0.18 

0.10 

0.04 

0.15 

E02IE0 

1.32 

0.67 

0.16 

0.14 

0.22 

0.19 

0.02 

EaxIEo 

3.67 

2.58 

1.16 

0.48 

0.17 

0.04 

0.01 

^Ai  ^Eo 

0.69 

0.25 

0.09 

0.05 

0.05 

0.04 

0.04 
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later  calculations  we  plan  to  do  only  a  quadreint  of  the  present  simiilations,  i.e.,  put  the 
expanding  debris  in  one  corner.  We  will  also  use  an  adjustable  grid  to  allow  the  shock 
to  move  out  of  the  system.  How  to  model  the  important  effect  of  the  recovery  of  the 
ambient  magnetic  field  in  this  case,  however,  remains  unclear  at  this  time.  Typically,  in 
these  calculations  we  use  200  x  200  cells  with  160000  background  amd  160000  debris  ions. 
In  this  section  we  show  results  of  three  different  calculations,  with  debris  Mach  number 
Md  =  VdI'^a  =  8,  0.2,  and  1,  respectively.  Scaling  relations  obtained  from  these  amd 
other  simulations  are  discussed  in  the  next  subsection. 

We  begin  with  the  Mq  =  8  case.  Representative  results  are  displayed  in  Figures  3.2- 
3.7.  Figure  3.2  shows  contour  plots  of  the  total  ion  density  at  various  times.  At  t  =  0 
we  can  see  only  the  debris  ions  (n£)o/no  =  32).  As  time  proceeds,  the  debris  density 
drops  and  the  outermost  contours  represent  air  ions  that  are  beginning  to  be  compressed 
by  the  expanding  debris.  By  Q.it  =  4.8,  the  outer  dark  edge  represents  the  front  of  an 
expanding  shock  in  the  air.  By  =  7.2  the  debris  ions,  whose  density  is  represented 
by  the  innermost,  wavy  contours,  are  essentizilly  stopped,  while  the  shock  continues  to 
expend  and  eventually  =  12)  reflects  from  the  ends  of  the  system.  Corresponding  3-D 
perspectives  of  the  magnetic  field  are  displayed  in  Figure  3.3  at  comparable  times.  The 
high  Mtch  number  expemsion  rapidly  compresses  the  magnetic  field,  which  then  expands 
out  as  a  collisionless  shock  wave.  The  region  behind  the  shock  remains  highly  ttirbulent. 
Again,  complex  interactions  occur  at  the  boundaries,  as  the  shock  arrives  and  reflects. 

Figxire  3.4  shows  a  contour  plot  of  the  magnetic  field  at  Qit  =  6.  The  smooth  outer 
contours  of  the  shock  region  and  the  more  structured  contours  in  the  region  of  the  debris 
are  evident.  Figure  3.5  displays  slices  in  x  of  the  magnetic  field  and  u*  —  i  phase  space  for 
both  the  air  and  debris  ions.  The  magnetic  field  profiles  show  the  expanding  shock  with  a 
large  overshoot,  sepeirating  out  from  the  driving  piston.  The  shock  has  evidently  energized 
the  air  ions,  by  the  reflection/gyration  process  discussed  in  the  previous  section.  By  this 
time  {Qit  =  6),  the  debris  has  stopped.  As  Table  3.1  shows,  the  equal  mass  radius  in  this 
case  is  12.8c/a;,-,  but  because  the  ion  gyroradius  is  so  large  {pi  =  Vjp/fl,  =  8),  the  debris 
actually  extends  out  as  far  as  Rm  -|-  2pi  ~  Z0c/u>i.  Most  of  the  debris  ions  have  stopped 
and  not  thermzdized  much.  Some  debris  ions,  however,  are  energized  diuring  the  debris/air 
interaction. 

In  Figure  3.6  are  displayed  3-D  perspective  plots  of  the  debris  density,  the  debris  energy 
[Ed  =  0.5mo(v^  +  Vy)],  and  the  air  ion  energy  at  later  time  =  16)  for  the  szune  case. 
All  three  quantities  are  chzuracterized  by  very  spiky  featxires  that  indicate  the  highly  tur¬ 
bulent  nature  of  the  interaction.  The  debris  ions,  however,  remain  fairly  well  contained  at 
about  the  same  radius  as  they  were  at  earlier  times.  Again,  this  is  because  the  present  cal¬ 
culations  do  not  model  well  the  recollapse  of  the  ambient  magnetic  field,  i.e.,  the  recollapse 
is  accompanied  by  the  reflected  shock  wave.  The  debris  energy  is  also  highly  structured, 
with  a  more  energetic  debris  core  and  a  more  highly  variable  energy  in  the  outer  halo  re¬ 
gion.  The  air  ion  energy  is  also  structured,  with  most  of  the  energization  occurring  in  the 
outer  regions,  as  expected.  In  order  to  show  these  effects  in  a  more  quantitative  fashion, 
narrow  slices  in  y  of  these  sjime  quantities  at  the  same  time  are  plotted  as  a  function  of 
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Uit  =  0 


na  =  7.2 


Figure  3.3-  3-D  perspectives  of  the  magnetic  field  at  the  same  time  for  the  same  run. 
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Figure  3.6.  3-D  perspectives  of  the  debris  density,  debris  energy,  and  air  ion  energy  at 
Qit  —  16  for  the  same  run. 


X  in  Figure  3.7.  The  spiky  nature  of  the  debris  ion  density  (normalized  in  terms  of  the 
peak  debris  density)  is  especially  evident  in  the  outer  regions.  Near  the  density 

remains  small;  better  collapse  of  the  magnetic  field  would  tend  cc  impress  the  density, 
giving  rise  to  a  more  expected  monotonic  falloff.  A  fraction  of  the  density  remains  '-i  a 
low  density  halo.  The  second  panel  in  Figure  3.7  shows  the  debris  energy,  normalized  to 
the  initial  debris  kinetic  energy.  As  with  the  density,  the  debris  ener^-  profile  can  be  sub¬ 
divided  into  an  iimer  enegetic  core  at  r  <  ri  and  an  outer  halo  of  less  energetic  debris  ions 
(ri  <  r  <  r2).  The  corresponding  air  ion  energy  is  plotted  in  the  bottom  panel.  The  inner 
region  is  essentially  devoid  of  ciir  ions.  Air  ions  outside  of  this  region  have  been  energized 
to  a  significant  fraction  of  the  initied  kinetic  energy  by  the  shock. 

In  contrcist,  we  plot  in  Figures  3.8-3.11  corresponding  quantities  from  a  low  Mach  number 
expansion  at  Mjrj  =  0.2.  Because  similar  low  Mach  number  cases  have  been  discussed 
previously  in  Section  1  in  relation  to  VHANE  bursts,  we  limit  the  discussion  here  to  a  few 
plots  to  compare  with  the  Md  =  8  case.  In  Figure  3.8  is  displayed  a  contour  plot  of  Bz  at 
flit  =  24.  (Note  that  this  calculation  is  done  in  a  64  x  64c/u;i  system  with  correspondingly 
fewer  cells  and  macroparticles.)  The  plot  shows  the  formation  of  a  magnetic  cavity  with 
weak  magnetosonic  waves  propagating  outwzird.  A  profile  of  the  magnetic  field  for  this  case 
is  shown  in  the  top  panel  of  Figure  3.9.  The  cavity  and  the  wezdc  magnetic  compression 
ahead  of  the  cavity  are  clearly  seen.  The  bottom  peinels  show  —  x  phase  space  for  the 
air  and  debris  ions  at  the  szune  time.  The  cavity  formation  has  led  to  some  slight  heating 
of  the  air  ions  at  the  edge  of  the  cavity.  The  debris  ions  zxe  still  expanding  at  this  time 
(Slit  =  24).  With  the  magnetic  confinement  radius  in  this  case,  Rb  =  12.8c/wi,  the  debris 
continues  to  expand  until  f2j<  ~  60. 

Three  dimensional  perspectives  of  the  debris  density,  debris  energy,  and  air  ion  energy  are 
shown  in  Figure  3.10  at  a  much  later  time,  =  160,  when  the  debris  heis  collapsed  with 
the  magnetic  field.  The  debris  remains  well  confined  in  this  case,  although  the  energy  of  the 
debris  is  quite  variable.  Unlike  the  high  Md  case,  here  the  air  ions  that  become  strongly 
energized  are  those  that  collapse  with  the  debris.  The  air  ions  outside  the  mmn  cloud, 
which  interact  with  weak  magnetosonic  disturbances,  are  not  heated  very  much.  These 
results  are  reenforced  in  Figure  3.11,  which  show  corresponding  cuts  of  these  distributions 
along  X.  The  debris  is  well  contained  within  a  radius  ri ;  a  small  fraction  resides  in  a  narrow 
hzJo  region,  ri  <  r  <  r2.  Air  ions  in  the  core  region  are  strongly  heated.  Less  air  heating 
occurs  in  the  halo  region,  even  less  so  at  larger  radii. 

Finally,  in  Figtires  3.12-3.17  we  show  results  from  an  intermediate  case,  with  Md  =  1- 
Figure  3.12,  like  Figures  3.4  and  3.8,  shows  a  2-D  contour  plot  of  Bz  at  fljt  =  32.  The  inner 
contours  correspond  to  the  location  of  the  diamagnetic  cavity,  while  the  outer  contour  is  the 
edge  of  a  weak  expanding  shock.  Figure  13  shows  the  corresponding  density  contours.  The 
inner  contours  represent  the  debris,  which  has  essentially  stopped  at  this  time,  compressed 
into  a  thin  shell,  and  structured.  The  outer  contours  are  £iir  ions  behind  the  shock  front. 
Figure  3.14  (again  like  Figures  3.5  and  3.9)  show  profiles  of  the  magnetic  field  (again  the 
cavity  and  shock  en-e  clearly  seen)  and  —  x  pheise  space  for  the  air  eind  debris  ions. 
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Figiire  3.7.  Slices  in  x  of  the  debris  density,  debris  energy,  and  air  ion  energy  at  the  same 
time  as  in  the  previous  figure. 


Figure  3.9.  Slices  in  x  of  and  Vg  —  x  phase  space  for  the  air  and 
time  as  in  the  previous  fig\ire. 
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Figure  3.11.  Slices  in  x  of  the  debris  density,  debris  energy,  and  air  ion  energy  at  the  same 
time  as  in  the  previous  figure. 
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Figure  3.13.  2-D  contour  plot  of  the  ion  density  at  the  same  time. 
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Figure  3.14.  Slices  in  x  of  and  Vs  —  x  phase  space  for  the  air  and  debris  ions 
tame  time  as  in  the  previous  figure. 


59 


The  air  ions  are  wezdcly  heated  behind  the  shock.  The  debris  ions  have  actually  begun  to 
collapse  (ui  <  0  at  the  right  edge  of  the  cloud)  with  the  magnetic  field  at  this  time. 

Three  dimensional  perspectives  and  radial  slices  of  the  debris  density,  debris  energy,  and 
air  energy  are  displayed  in  Figures  3.15-3.16.  Like  Figures  3.10-3.11  for  the  Md  =  0.2 
case,  the  debris  remains  well  confined  and  the  air  ions  trapped  in  the  cloud  axe  strongly 
energized.  Like  Figures  3. 6-3. 7  for  the  Md  =  8  case,  air  ions  outside  the  main  debris  cloud 
have  been  heated  by  the  shock.  And  overall,  there  are  large  fluctuations  of  each  quantity 
in  the  cloud. 

In  Figure  17  are  plotted  energy  distributions  for  the  debris  ions  for  the  same  Md  =  1 
case  at  various  times.  Initially,  the  debris  energy  peaks  around  Eg  with  some  thermal 
spread.  As  the  debris  expands  and  stops,  the  spectrum  becomes  peaked  at  F  ~  0  auid 
is  fairly  cold.  As  the  magnetic  field  recovers  (Qit  —  48),  the  debris  ions  au-e  heated. 
Eventually  =  120),  the  spectrum  consists  of  dense,  moderately  hot  component  and  a 
more  energetic  tail,  extending  out  to  F  ~  3Fo. 

Figure  3.18  displays  debris  ions  at  the  end  of  the  simulations  for  the  three  cases  we  have 
disciissed  in  detail:  Md  =  0.2,  1,  amd  8  as  well  as  for  Md  =  2.  The  Md  =  0.2  amd  1  cases 
show  the  debris  has  collapsed  to  a  relatively  small,  highly  structured  cloud.  The  Md  =  2 
case  also  exhibits  structure,  but  the  debris  ions  lie  in  a  ring  distribution.  If  run  longer,  the 
debris  will  slowly  collapse  further.  In  the  Md  =  8  case  a  smadl  central  component  of  the 
debris  has  collapsed,  but  much  of  the  debris  remains  out  at  its  point  of  maximum  extent 
and  is  highly  structured.  In  this  case,  however,  the  ambient  field  does  not  recover  well,  so 
that  the  final  debris  size  is  not  resolved.  The  larger  spatial  extent  of  the  debris  in  this  case 
is  probably  also  due  to  the  large  ion  gyroradius  {pi  =  Moc/wj). 

3.3  Scaling  Plelations 

In  addition  to  the  three  cases  discussed  above,  we  have  carried  out  additioned  simulations 
of  the  same  type,  varying  only  the  debris  Mach  number.  The  results  of  these  simulations 
zure  summarized  in  Table  3.1.  For  each  case  we  give  the  initizd  debris  Mach  number, 
debris  density  {riDo/no),  debris  energy  (Fo/0.5moU^),  the  magnetic  confinement  radius 
Rb,  and  equzd  mass  radius  Rm.  At  the  end  of  the  run,  when  the  debris  has  collapsed, 
we  calculate,  as  in  Figures  3.7,  3.11,  and  3.16,  a  cloud  radius  rj  and  a  halo  radius  r2, 
normalized  to  the  smaller  oi  Rm  or  Rb,  the  debris  and  air  density  zind  energy  inside  the 
cloud  (rtDi,  Fdi,  Eai)  and  in  the  halo  (nD2»  Ed2^  normalized  to  noo, 

n„  and  Eo-  As  in  Section  1  of  this  report,  rj  is  defined  to  be  the  radial  position  where 
the  density  drops  to  20%  of  its  peak  value;  r2  is  the  maxim\im  radial  extent  of  the  debris. 
Again,  as  in  Figure  3.7  and  3.18  indicate,  at  higher  Md  the  debris  does  not  fully  collapse, 
so  the  final  size  of  the  debris  cloud  may  be  overestimated  zmd  the  energy  in  the  debris  may 
be  underestimated. 

Nevertheless,  some  interesting  conclusions  can  be  drawn  from  the  results  in  the  table, 
which  have  been  plotted  in  Figmes  3.19-3.20.  The  debris  radii,  rj  and  r2  normalized  by 
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Figtire  3.19.  Results  of  simulations  with  various  Mp  showing:  (top  panel)  normalized 
radii,  rj  and  rj  (normalized  to  r,  =  confinement  radius  =  smaller  of  Rg  and  Rm); 
(middle  panel)  debris  density,  npi,  np^;  and  (bottom  panel)  air  density,  tiai,  ^A2- 
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Figure  3.20.  Results  of  simulations  with  various  Md  showing:  (top  panel)  debris  energy, 
Ed2i  Mid  (bottom  panel)  air  energy,  E^. 
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r,,  the  smaller  of  the  two  stopping  radii,  are  plotted  versus  M©  in  the  top  panel  of  Figure 
3.19.  Both  normalized  radii  are  essentially  constant  out  to  Md  —  1-  For  larger  Md,  the 
normalized  r2  remains  about  constant  (except  &t  Md  =  8).  This  one  larger  value  of  r2  and 
the  larger  values  of  rj  for  Mq  >  1  is  due  in  part  to  the  lack  of  the  magnetic  field  collapse 
eluded  to  earlier  and  to  the  larger  ion  gyroradii  at  larger  Md-  On  the  other  hand,  it  is 
also  possible  that  the  debris  confinement  radius  does  increase  at  higher  Md  just  because 
the  magnetic  field  plays  a  less  significant  role.  The  middle  panel  of  Figure  3.19  shows  the 
debris  density  tidi  inside  the  cloud  and  in  the  outer  haJo  (n£)2),  normalized  to  the  peak 
debris  density.  nD2  is  about  constant  out  to  Md  —  1,  then  decreases  at  larger  Md-  The 
faiioff  of  UDi  is  due  in  part  to  the  fact  that  the  initial  density  relative  to  the  bax;kground 
density  decreases  with  Md  (to  Md  =  1)  (Table  3.1].  At  higher  Md,  when  the  initiail 
peak  debris  density  is  constant,  riDi  decreases  because  ri  increases.  The  corresponding  air 
densities  are  shown  in  the  bottom  panel  of  Figure  3.19.  The  air  density  inside  the  debris 
cloud  is  roughly  constant  at  small  Md  and  decreases  at  large  Md,  because  the  expanding 
shock  at  high  Mach  numbers  tends  to  expel  the  air  ions.  The  air  density  outside  the  main 
debris  cloud  is  essentially  constant.  The  next  two  graphs  (Figure  3.20)  show  the  energy 
of  the  debris  (top  panel)  and  air  (bottom  panel)  ions  in  the  main  debris  cloud  and  in 
the  halo  region.  The  energy  tends  to  fall  off  monotonically  with  Md-  At  small  Md  the 
energy  can  exceed  the  initial  energy;  ions  left  in  the  cloud  are  heated  at  the  magnetic  field 
recompresses.  The  lower  energy  (~  lQ%Eo)  at  higher  Md  (>  1)  reflects  the  fact  that  much 
of  the  initial  energy  goes  into  heating  of  the  background  ions  outside  of  the  debris  cloud 
eind  in  making  large  amplitude  magnetic  disturbances. 

3.4  Summary 

In  this  section  we  have  discussed  debris  and  air  dyneimics  related  to  VHANE  and  HANE 
bursts.  In  particular,  we  have  examined  a  series  of  hybrid  simulations  in  which  the  debris 
Mach  number  is  varied  systematically.  We  showed  results  for  three  typical  cEises:  a  high 
Mach  number  expansion  that  corresponds  to  a  weak  HANE  burst,  a  low  M£w:h  number 
expansion  characteristic  of  a  VHANE  burst,  and  an  intermediate  case  at  Md  =  1-  We 
also  carried  out  a  number  of  other  rims  and  characterized  the  debris  and  air  properties  at 
the  end  of  the  calculation,  when  the  debris  and  magnetic  field  recollapsed,  as  a  function 
of  Md- 

Generally,  we  find  in  the  low  Md  expansion,  similar  to  those  discussed  in  Section  1  of  this 
report,  that  the  debris  remains  well  contained,  although  there  is  the  usual  structuring  of 
the  debris  when  it  is  stopped  by  the  magnetic  field.  On  recovery  of  the  magnetic  field 
the  debris  also  collapses,  maintaining  its  highly  structured  nature,  and  is  heated  in  the 
process.  Air  ions  trapped  inside  the  debris  are  also  heated.  Outside  the  debris  cloud  there 
is  little  wave  activity  or  air  heating.  In  contrast,  at  higher  Md  a  shock  forms  as  the  debris 
expands,  energizing  the  air  ions  in  the  outer  region.  Constrmnts  of  the  present  calculations 
do  not  allow  the  correct  recollapse  of  the  magnetic  field,  so  the  debris  tends  to  remain  out 
to  a  radius  corresponding  to  the  equal  mass  radius.  Heating  of  the  debris  is  smaller  in  this 
case,  as  much  of  the  energy  goes  into  the  magnetic  field  disturbemces  and  heating  of  the 
air  ions.  The  intermediate  Md  =  1  ceise  shows  aspects  of  both  types  of  expzmsions.  Like 
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the  lower  Md  cases,  the  debris  structures  at  r  ~  Rb  —  Rmi  recollapses,  and  heats.  Like 
the  higher  Md  runs,  a  (weak)  expanding  shock  is  formed  and  heats  the  air  ions  to  some 
degree.  It  is  not  evident  here  that  Md  =  1  is  a  special  case  or  gives  rise  to  particularly 
large  magnetosonic  waves. 

In  this  study  we  have  begim  to  develop  scaling  relations  for  HANE/VHANE  bursts.  It 
appears  that  the  debris  stays  well  confined,  t.e.,  no  large  jetting  is  seen,  at  all  Mach 
numbers.  At  lower  Md,  the  final  radius  is  about  half  of  the  confinement  radius;  at  larger 
Md  the  final  radius  is  roughly  equal  to  the  confinement  radius.  It  has  also  been  shown 
that  the  energy  in  the  debris  and  air  ions  decreases  with  Md-  The  energy  remedning  in 
the  debris  decreases  to  about  10%  of  the  initial  debris  kinetic  energy  for  Md  ^  1,  as 
more  energy  goes  into  heating  of  the  baekgroimd  ions  and  in  the  generation  of  the  shock. 
Because  the  present  calculations  are  done  in  the  two  dimensions  perpendicvdar  to  the 
eimbient  magnetic  field  direction,  however,  we  cannot  say  how  much  debris  remains  in  the 
burst  region  when  dynaimics  along  B  are  included,  and  how  that  fraction  varies  with  Md- 
Later  calculations  will  include  such  paxcdlel  eflFects  and  allow  us  to  estimate  this  fraction. 
We  also  eventually  plan  to  compare  the  results  here  with  the  more  general  2-D  calculations 
of  Brecht  [unpublished,  1991]  and  Simonson  emd  Hewett  [1991]. 
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4.  Parallel  Streaming  Deposition  Model 

4.1  Introduction 

So  far,  we  have  discussed  the  case  of  debris-air  dynamics  perpendicular  to  the  magnetic 
field  in  both  the  HANE  and  VHANE  regime.  When  debris  streams  out  perpendicular  to 
B,  it  is  stopped  by  air  ions  in  the  HANE  regime,  and  the  magnetic  field  in  the  VHANE 
case.  Some  debris  ions  are  deposited  in  the  ionosphere,  either  because  they  stream  down 
(or  up)  field  lines  or  because  they  are  turned  from  streaming  perpendicular  to  the  field 
to  the  paredlel  direction,  as  discussed  in  Section  2.  Many  more  air  ions  are  energized  and 
deposited  in  the  ionosphere  as  well,  creating  the  kinetic  energy  patches.  How  far  the  ions 
go  into  the  ionosphere  and  where  they  deposit  their  energy  is  generally  treated  as  a  simple 
single  particle  collisional  process.  At  lower  altitudes  the  background  air  density  increases 
and  the  mean  free  path  of  the  ions  thus  becomes  shorter.  The  altitudes  at  which  such 
collisional  processes  become  important  are  typically  below  200  km;  above  this  height  the 
mean  free  path  is  so  large  that  the  ions  essentially  free  stream.  For  example,  Kilb  [1990] 
shows  that  an  Fe"*"  ion  with  a  velocity  of  270  km/sec  has  a  meam  free  path  of  38000  km 
when  interacting  with  the  atmosphere  at  400  km. 

However,  it  is  well  known  that  anomalous  wave-particle  processes  in  a  plasma  can  reduce 
this  streaming  distance  significantly.  For  example,  ions  backstreaming  from  the  Earth’s 
bow  shock  into  the  solar  wind  have  a  classical  mean  free  path  of  about  10®  km.  In  situ 
satellite  measurements  indicate  that  such  ions  are  actually  stopped  over  a  distance  on  the 
order  of  10®  km  [Bonifaizi  and  Moreno,  1981].  This  anomalously  short  stopping  distance 
occurs  because  the  backstreaming  ions  generate  electromagnetic  instabilities  which  produce 
waves  that  in  turn  scatter  the  ions  [see  Gary  (1991)  for  a  recent  review).  These  same  kinds 
of  electromagnetic  instabilities  have  been  observed  in  hybrid  simulations  of  HANEs  [Brecht 
and  Thomas,  1988].  The  issue  to  be  addressed  here  is  how  one  can  estimate  the  effective 
stopping  distance  due  to  such  instabilities  and  whether  such  am  effect  can  be  important  for 
HANE  and  VHANE  bursts.  In  Section  4.2  we  review  the  basic  physics  and  derive  some 
elementary  scaling  relations.  Then  in  Section  4.3  we  apply  these  conditions  to  some  simple 
models  for  HANE  bursts.  Finally,  in  Section  4.4  we  draw  some  conclusions  and  suggest 
future  refinements  of  the  basic  model. 

4.2  Background 

Electromagnetic  instabilities  resulting  from  the  relative  strezuning  of  two  ion  components 
have  been  studied  for  a  long  time  [Gary  et  ai,  1984;  Winske  and  Leroy,  1984;  Gary  et  ai, 
1986;  Gary,  1991].  Generally,  these  studies  have  involved  homogeneous  beams  in  space. 
The  linezu-  [Gary  ei  al.,  1984]  and  quasi-linear  [Winske  and  Leroy,  1984]  properties  of  the 
instabilities  have  been  studied  analytically  and  numerically,  while  the  nonlinear  behavior 
have  been  deduced  from  hybrid  simulation  studies  [Winske  and  Leroy,  1984;  Gary  et  ai, 
1986;  Winske  and  Quest,  1986;  Thomas  emd  Brecht,  1986].  More  recently,  the  instability 
has  been  studied  in  time  [Miller  et  al,  1991]  and  space  [Onsager  et  al.^  1991]  varying 
situations.  These  electromagnetic  instabilities  are  particularly  interesting  in  the  HANE 


69 


regime  because  they  are  quite  robust,  i.e.,  they  do  no  depend  strongly  on  the  details  of  the 
velocity  distribution  aind  they  persist  in  a  highly  turbulent  regime,  when  even  the  direction 
of  the  magnetic  field  is  not  well  defined. 

We  consider  a  field  aligned  ion  beam,  density  and  ion  mass  mj,  streaming  with  ve¬ 
locity  Vfc  relative  to  a  background  plasma,  ion  density  and  mass  rrij,  along  an  am¬ 
bient  magnetic  field  Bo.  We  define  the  total  density  rio  =  ni,  +  Ui  and  the  Alfven 
speed  va  =  Bo/i^nnomiY^^ .  We  generally  consider  the  plasma  components  as  cold, 
i.e.,  =  SnrioTj/B'^  <  1  where  j  is  the  plasma  species  (beam,  background,  or  electrons). 

In  this  situation  two  electromagnetic  ion  beam  instabilities  can  occur.  One  instability 
involves  the  (anomalous)  cyclotron  resonance  of  the  Doppler  shifted  wave  with  the  bezun, 
i.e., 

‘i'r  -  fcVJ,  ~  -fifc  (4.1) 

where  Clj  =  eBo/mjC  is  the  cyclotron  frequency  of  the  j-th  species.  This  instability  will 
be  referred  to  as  the  “resonant”  mode.  The  other  instability  is  a  nonresonant  (i.e.,  fluid) 
interaction  that  involves  the  presence  of  the  beam  to  give  an  effective  T||/T_l  >  1.  This 
mode,  termed  here  tne  “noru'esonant  instability”  is  thus  a  fire-hose-like  instability.  There 
is  no  corresponding  resonance  condition  like  Eq.  (4.1)  for  the  nonresonant  instability. 

The  two  instabilities  have  distinct  Unear  properties,  which  are  exhibited  in  Figures  4.1- 
4.4.  Figure  4.1  makes  the  point  discussed  above:  namely,  that  the  beam  ions  are  in 
cyclotron  resonance  with  the  waves  for  the  resonance  mode,  while  waves  generated  by 
the  nonresonant  instabiUty  do  not  fall  in  velocity  space  where  either  ion  species  resides. 
Generally,  the  phase  velocity  of  the  resonant  waves  Ues  in  the  vicinity  of  the  beam  ions, 
while  that  of  the  nonresonamt  instability  lies  close  to  the  backgroimd  ions.  Figure  4.2 
shows  growth  rates  for  the  resonant  (R)  and  nonresonant  (NR)  instabiUties  maximized  over 
wavenumber  versus  b'  im  density  nj/no  for  VhfvA  =  10  and  =  mj.  The  nonresonant 
instability  has  a  finit :  threshold  for  growth.  Above  threshold  the  growth  rate  increases 
linearly  with  ruing.  On  the  other  hand,  the  resonant  instability  has  no  threshold:  positive 
growth  persists  to  Uf,  —*  0.  However,  the  maximum  growth  rate  increases  more  slowly 
with  beam  density  so  that  eventually  {n(,/no  >  0.08)  the  growth  rate  of  the  nonresonant 
instability  is  Izirger.  Figure  4.3  shows  maximum  growth  rates  eis  a  function  of  for 

nh/uo  =  0.1.  Both  instabilities  has  a  finite  threshold,  with  that  of  the  resonant  mode 
being  smaller.  Above  threshold  7  increases  Unetirly  with  Vj/v,4  for  the  nonresonzmt  mode 
to  Vi/VA  15;  thereafter  it  increeises  at  a  slower  rate.  On  the  other  hand,  the  growth  rate 
for  the  resonant  mode  increases  very  fast  above  threshold  to  a  maximum  value;  at  higher 
Vhl^A  the  maximum  growth  rate  is  essentially  constant.  Figure  4.4  presents  threshold 
values  for  the  instabilities  as  a  fimction  of  Vh/vA  and  ni/ng.  Generally,  the  resonemt 
instability  (solid  curve)  has  a  lower  threshold,  Vt/vA  <  2,  which  does  not  vary  much  with 
beam  density.  (The  dotted  Une  corresponds  to  conditions  when  the  two  modes  have  equal 
growth  rates;  above  the  line  the  noruresonant  mode  is  more  unstable.) 
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RIGHT  HAND  NONRESONANT 


Figtire  4.1.  Schematic  of  background  and  beam  velocity  distributions  along  with  the  loca¬ 
tion  of  cyclotron  resonances  for  the  two  instabilities. 
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rates  for  the  resonant  (R)  and  nonresonant  (NR)  instabilities  maxi- 
nvunber  versus  n»/no  for  V»/vyi  =  10- 
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Figure  4.3.  Growth  rates  maximized  over  wavenumber  for  the  two  instabilities  as  a  function 
of  Vhfvji  for  riifno  —  0.1. 
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rih/no 


Figure  4.4.  Thresholds  for  the  resonant  (solid  curve)  and  nonresonant  (dashed  curve) 
instabilities  as  a  function  of  Vj/oa  and  nj/rio.  The  dotted  curve  is  where  the  two 
instabilities  has  equal  growth  rates;  above  the  curve  the  nonresonant  instability  dom¬ 
inates. 
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We  can  summarize  the  linear  properties  of  the  two  instabilities  as  follows.  The  growth 
rates  of  the  two  modes  vary  as: 


(4.2) 

2''no  ''mj,'' 

7  ^  1  nj  Vj 

Cli  ~  2  rio  va 

(4.3) 

s  can  be  expressed  as: 

Vfc  ~  1  —  2va 

(4.4) 

Vb  > 

Tib 

(4.5) 

Again,  this  result  follows  from  Figure  4.4. 


In  addition  to  these  linear  properties  the  nonlinear  evolution  of  these  two  instabilities  has 
been  studied  extensively  via  hybrid  simulations  in  both  1-D  [Winske  aind  Leroy,  1984;  Gary 
et  ai,  1986]  and  in  2-D  [Winske  and  Quest,  1986;  Thomas  jind  Brecht,  1986].  Figure  4.5 
shows  an  example  of  a  1-D  calculation.  Here  Vb/vA  is  teiken  to  be  10  and  nb/uo  =  0.02, 
i.e.,  in  this  case  the  resonant  instability  dominates.  The  left  panels  of  the  figure  show 
Vx  —  X  phase  space  for  the  two  ion  species  at  vauious  times.  The  right  panels  show  profiles 
in  X  of  one  component  of  the  magnetic  field  at  the  same  times.  Initially,  there  are  two  well 
defined  ion  streams  and  no  fluctuations.  As  the  instability  is  excited,  waves  grow  and  the 
beams  begin  to  couple  together.  Times  are  expressed  here  in  terms  of  the  inverse  growth 
time.  By  t  =  5/7  coupling  of  the  two  beams  has  begun,  and  by  <  ~  IO/7  the  beaims  have 
completely  merged. 


One  issue  that  is  particularly  important  for  the  HANE  case  is  the  effect  of  a  finite  length 
beam  (as  we  will  see  later),  which  occurs  as  the  debris  ions  compress  into  a  thin  coupling 
shell.  Figure  4.6  shows  a  simulation  of  a  finite  length  beam  segment  interacting  with  a 
uniform  background  plasma  [Onsager  et  al.,  1991].  One  sees  that  very  quickly  (fl.t  =  10) 
there  is  already  the  beginning  of  significant  slowing  of  the  beeim  and  background  pleismas 
and  eventually  both  are  strongly  coupled.  In  this  case  rib/rio  =  0.4  emd  Vb/vA  =  10,  so 
that  Eq.  (4.2)  imphes  'y/Cli  —  0.4  and  Qit  =  10  corresponds  to  jt  =  4.  Detailed  analysis  of 
such  simulations  indicate  that  while  the  nonresonant  instability  has  a  larger  hnezir  growth 
rate,  the  resonant  mode  is  actually  more  important  because  the  resonant  mode  rem£uns  in 
contact  with  the  finite  length  beam  for  a  longer  period  of  time.  Based  on  these  simulation 
results  and  linear  theory,  we  can  then  construct  a  coupling  model. 


4.3  HANE  Coupling  Model 


The  HANE  anomalous  debris  coupling  model  is  based  on  Eqs.  4. 2-4. 5  and  the  above 
discussion  concerning  the  simulations.  The  basic  assumption  of  the  model  is  that  when  the 
debris  density  exceeds  the  ambient  air  density,  the  relevant  instability  is  the  nonresonant 
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Figure  4.6.  Results  of  a  1-D  hybrid  simulation  with  a  finite  length  ion  beam  (n^/n,,  = 
0.4,  Vb/vA  =  10)  showing  rapid  coupling  of  the  two  beams  due  to  the  resonant  insta¬ 
bility. 
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mode.  In  this  case  the  debris  acts  like  the  backgroimd  plasma  and  the  nonresonant  mode 
is  most  important  because  it  propagates  with  a  phase  velocity  near  that  of  the  debris.  The 
resonant  instability  is  more  relevant  when  either  the  nonresonant  instability  is  stable  (it 
hzis  a  higher  threshold  than  the  resonant  mode)  or  when  the  debris  density  is  less  them 
that  of  the  ambient  air.  In  either  case  we  compute  the  (local)  growth  rate  7  and  take  the 
coupling  time  to  be  5/7.  In  terms  of  equations,  the  model  is  thus: 


1.  When  no  >  and 

Vo  > 

riA  no 

the  nonresonant  mode  dominates.  Its  growth  rate  is 


7 


1/2  1/2 
’'A  "g 

riA  +  no  vy^  mo 


where  va  is  computed  using  the  backgroimd  air  density  and  mass. 


(4.6) 


(4.7) 


2.  When  no  <n a  and  Vo  >  1-5va,  the  resonant  mode  dominates  with 


1.  np 

2  ua  +  np 


mp 


(4.8) 


In  both  cases  the  coupling  time  is 

Tc  ~  5/7 

and  the  anomalous  coupling  distance  is  roughly 

Xan  ~  5Vb/7 


(4.9) 


(4.10) 


This  simple  model  does  not  necessarily  give  a  smooth  transition  from  one  regime  to  another, 
nor  does  it  attempt  to  work  out  the  continuous  slowing  of  the  debris.  As  we  show  next, 
we  axe  only  interested  at  this  point  in  comparing  these  coupling  times  and  distances  with 
their  collisional  counterparts  to  look  for  large  differences  between  the  two  rates. 

We  assume  a  1  MT  burst  yielding  10*®  J  of  kinetic  energy,  i.e.,  500  kg  of  mass  with  velocity 
Vp  =  2000  km/sec.  Taking  the  debris  ion  atomic  mass  A=30  then  yields  Np  =  10^®  debris 
ions.  We  shall  assume  two  spatial  models  for  the  debris  that  represent  extreme  limits.  First, 
we  assume  that  the  debris  expzinds  into  a  sphere  of  radius  R  with  uniform  density,  i.e.. 

No  =  yR’nD(i2)  (4.11) 

or 

o  4  X 

noiR)  =  -  V 

/t'*(km) 
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We  consider  a  burst  at  800  km,  roughly  halfway  between  STARFISH  altitudes  (~  400  km) 
and  the  VHANE  regime  (~  1200  km).  We  take  air  mass  densities  from  Kilb  [1990]  along 
with  collisional  coupling  lengths,  X^.  Taking  the  magnetic  field  to  be  0.3G  and  the  air 
atomic  mass  to  be  14,  we  find  that  the  ambient  air  number  density  is 

riA  =  4.3  X  lO^^pfc  cm~^  (4.13) 


where  pb  is  the  air  mass  density  in  gm  per  cm  The  Alfven  speed  is 

va  ~  1.75  X  lO^n^^'^^  km/sec 
At  altitude  h,  the  debris  density  is  thus 


noih) 


2.4  X  10^2 

(800  -  hy 


(4.14) 


(4.15) 


In  Table  4.1  are  tabulated  values  for  P6,  n^t,  n/j,  and  Va  at  various  altitudes  h.  We  see  that 
in  this  model  the  debris  density  is  always  less  than  the  aiir  density.  In  this  case  we  use  Eq. 
(4.8)  to  evaluate  the  growth  rate  and  (4.10)  to  estimate  the  anomalous  coupling  distance. 
The  last  column  in  the  table  is  the  collisioral  coupling  distance.  As  one  might  expect  from 
this  model,  at  higher  altitudes  the  debris  density  is  Izurger,  giving  rise  to  larger  growth  rates 
and  hence  shorter  coupling  distances.  Above  ~  330  km  the  anomalous  coupling  distance 
is  shorter  than  the  collisional  stopping  length.  At  600  km,  where  the  collisional  mean  free 
path  is  essentially  infinite,  the  anomalous  stopping  disteince  is  about  400  km.  Because  the 
anomalous  stopping  distance  increases  at  lower  altitude,  Xan  ~  400  km  is  a  lower  limit. 
From  Table  4.1  we  see  that  the  collisional  stopping  distance  at  200  km  (i.e.,  600  -  400  km) 
is  also  small.  Thus,  one  would  expect  that  the  overall  collisional  eind  anomalous  stopping 
distances  are  about  the  seune  in  this  case.  However,  the  detzuls  of  how  the  ions  slow  down 
amd  deposit  their  energy  along  the  trajectory  may  be  different. 


The  second  limiting  Ccise  we  consider  assumes  that  the  debris  compresses  into  a  thin  cou¬ 
pling  shell  of  thickness  p,  =  Vb/flo  ~  20  km.  In  this  case  the  debris  density  falls  off  at  a 
much  slower  rate 

Nd  =  4nR^nD{R)A  (4.161 

with  A  ~  Pi  so  that  for  an  800  km  burst 


no 


4  X  10*® 

(800  -  hy 


cm 


-3 


(4.17) 


As  we  show  in  Table  4.2,  the  debris  density  now  exceeds  the  air  density  at  all  altitudes 
down  to  200  km.  In  this  case  the  nonresonant  instability  [Eq.  (4.7)]  dominates.  At  lower 
altitudes  the  debris  density  drops,  while  the  air  density  increases,  which  increases  the 
ratio  of  n^/(n^  4-  np).  The  ratio  of  Vq/va  also  rises  rapidly,  so  that  the  growth  rate 
increases.  Consequently,  the  anomalous  stopping  distemce  becomes  shorter.  For  ah  .alues 
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Table  4.1.  Comp^l^ison  of  anomalous  and  collisional  stopping  distances  for  a  burst  at  800 
’an  with  a  uniform  debris  distribution  [Eq.  (4.11)]. 


h{km) 

Pbig/cm^) 

nj\(cm 

nDicm~^) 

VAikm/s) 

7(sec  ’) 

Xa„(km) 

^c/(km) 

200 

2.5  X  10-^3 

1.1  X  10’° 

1.1  X  10^ 

1.7 

0.6 

16100 

55 

250 

6.3  X  10-^“* 

2.7  X  10° 

1.4  X  lO'* 

3.4 

1.1 

9300 

320 

300 

2.1  X  10-^^ 

9.0  X  10» 

1.9  X  10^ 

5.8 

1.7 

5900 

1050 

400 

9.2  X  10-^® 

4.0  X  10^ 

3.8  X  10^ 

28 

6.0 

1660 

4  X  10^ 

600 

1.0  X  10-’® 

4.3  X  10° 

3.0  X  10^ 

84 

25. 

400 

oc 
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Table  4.2.  Compaxison  of  anomalous  and  collisional  stopping  distances  for  a  burst  at  800 
km  with  a  shell  debris  distribution  [Eq.  (4.16)]. 


h{kTn) 

Pbig/cm^) 

riAicm 

n£,(cm-°) 
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55 
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50 
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of  h,  Xan  is  less  them  Xei  by  about  three  orders  of  magnitude.  Thus,  even  if  we  double  the 
coupling  times  [Eq.  (4.10)]  and  increase  the  thickness  of  the  debris  shell  by  a  factor  of  5, 
Xan  will  still  be  less  than  Xd  by  two  orders  of  magitude.  We  should  also  remark  that  the 
analytic  approximations  for  the  growth  rates  may  not  be  valid  for  very  extreme  parameters 
{VofvA  ^1)  such  that  the  coupling  distance  becomes  compzirable  to  the  wavelength  of 
the  mode. 

4.4  Summary 

In  this  section  we  have  used  well  known  properties  of  electromagnetic  ion  beam  instabilities 
to  derive  a  simple  model  for  anomalous  debris-air  coupling  when  the  debris  streams  along 
field  fines.  The  model  is  based  on  the  local  debris  and  air  densities  and  the  relative 
streaming  velocity  compared  to  the  local  Alfven  speed.  When  the  debris  density  exceeds 
that  of  the  air,  the  nonresonant  instability  dominates,  with  a  threshold  condition  for  growth 
given  by  Eq.  (4.6)  and  a  growth  rate  by  Eq.  (4.7).  This  type  of  situation  might  occur  in 
a  HANE  burst  if  the  debris  velocity  spread  is  fairly  smadl  and  the  debris  streams  out  in  a 
narrow  shell.  In  this  case  the  anomalous  coupling  distance  decreases  at  lower  altitudes,  as 
does  the  usual  collisional  stopping  distance.  These  two  effects  then  can  work  together  to 
deposit  the  energy  at  higher  altitudes. 

A  more  likely  situation,  however,  is  that  the  debris  will  spread  out  over  a  larger  spatial 
region  with  a  lower  density,  less  than  that  of  the  air  (except  for  very  high  altitude  bursts). 
In  this  case  the  resonant  electromagnetic  ion  beam  instability  provides  the  coupling.  As 
the  ratio  of  debris  to  air  density  decreases  at  lower  altitudes,  this  type  of  coupling  becomes 
less  effective,  unlike  collisional  processes. 

Thus  fau",  we  have  shown  that  under  some  conditions  anomalous  slowing  of  energetic  debris 
as  well  as  of  air  ions  can  yield  stopping  dist2inces  (mean  free  paths)  much  shorter  than 
those  obtained  by  collisional  effects.  This  has  been  demonstrated  for  specific  point  by 
point  comparisons.  What  is  actually  needed,  however,  is  a  formula  valid  for  the  entire 
path  length  of  an  energetic  ion  to  give  the  integrated,  rather  than  the  local,  effect.  When 
the  slowing  down  rate  are  small,  one  could  imagine  using  a  WKB-fike  approach  to  estimate 
the  cumulative  effect  of  the  slowing  of  an  energetic  ion.  And  one  may  wish  to  consider 
the  ions  as  deposited  according  to  where  they  loose  their  energy,  rather  thain  to  their  final 
stopping  point  [Sappenfield,  private  communication,  1990).  To  investigate  these  processes 
more  accurately,  one  should  carry  out  simulations  of  debris  ions  with  a  finite  spatial  spread 
interacting  with  an  air  background  that  contains  an  exponentieJly  increasing  density.  Such 
czJculations  would  then  allow  comparison  with  the  simple  formulais  derived  here  as  well  as 
testing  the  sensitivity  of  parameters  of  the  debris  deposition  model.  We  defer  this  more 
detailed  study  to  the  future. 
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5.  Multibursts 


5.1  Introduction 

The  question  of  nuclear  phenomenology  related  to  m\iltiple  nuclear  bursts  at  high  altitudes 
remains  largely  unknown.  The  systems  codes  generally  treat  the  effects  of  each  burst  sep¬ 
arately,  but  include  atmospheric  heave  due  to  earlier  bursts  [Sappenfield,  1987].  However, 
the  dynamics  of  debris  and  air  interactions  in  such  situations  has  not  been  studied  in  great 
detail.  In  this  section  we  present  a  few  preliminziry  results  of  kinetic  simulations  of  two 
simultaneous  bursts  in  the  VHANE  regime.  This  regime  is  chosen  because  the  dynamics 
are  simpler  and  better  tmderstood.  In  particular,  we  will  show  jin  example  of  a  full  par¬ 
ticle  calculation  that  indicates  both  electron  and  ion  heating  exceeding  that  from  a  single 
burst.  We  also  present  an  example  of  a  hybrid  code  run  with  a  backgrotmd  plasma  that 
shows  different  ion  dynamics.  We  feel  that  multibursts  in  the  VHANE  regime  edso  offer  the 
possibility  of  doing  unique  ion/ion  interactions  on  the  AGNES  laser  facility,  as  discussed 
later. 

5.2  Particle  Simulation  of  a  Multiburst 

In  the  past  few  years  we  have  carried  out  particle  simulations  of  VHANEs  over  a  wide 
range  of  conditions  [Winske,  1987,  1988,  1990;  Winske  and  Jones,  1990).  Most  of  these 
calculations  have  been  done  in  the  large  gyroradius  limit  (pi/Rs  =ion  gyroradius/magnetic 
confinement  radius  >1)  and  show  the  development  of  a  structuring  instability  on  the 
surface  of  the  expanding  plasma  cloud  and  related  ion  and  electron  heating.  Here,  we 
carry  out  a  comparable  particle  simulation  of  two  simultaneous  binrsts.  In  this  case  we 
take  two  equal  bursts,  each  with  pi  =  Rb  =  Sc/wj,  expanding  in  2-D  perpendicular  to  an 
ambient  magnetic  field.  The  separation  between  the  two  bursts  in  this  case  is  also  Rb-  The 
results  of  the  calculation  are  shown  in  Figures  5. 1-5.3.  Figure  5.1  shows  the  debris  ions  in 
x-y  space  at  several  times.  The  ions  for  each  burst  are  displayed  in  separate  (left-right) 
panels.  As  in  a  single  burst,  the  ions  expand,  compress  into  a  shell,  amd  structure.  At  later 
times  the  inner  most  ions  contract,  while  the  ions  on  the  outer  edge  continue  to  expamd, 
leading  to  a  cloud  that  extends  out  beyond  Rb-  In  addition  to  these  single  burst  effects, 
there  is  obviously  an  interaction  going  on  where  the  two  clouds  intersect.  The  ions  in  this 
region  are  more  diffuse,  i.e.,  the  flutes  are  somewhat  more  elongated  and  the  ions  are  more 
scattered.  This  asymmetry  persists  to  late  times.  Note  in  particular  that  the  debris  from 
the  two  clouds  tends  to  mix  together. 

Figure  5.2  shows  corresponding  plots  for  the  electrons  from  the  two  bursts.  As  with 
the  ions,  one  sees  expansion,  compression,  and  structure.  However,  the  electrons  do  not 
mix  together.  This  is  not  unexpected,  as  the  electron  gyroradius  is  much  smaller.  Charge 
neutradity  is  maintained  as  the  electrons  associated  with  the  debris  ions  from  one  burst  just 
pair  up  with  other  debris  ions  at  late  times  as  the  ions  (with  larger  gyroradius)  intersect. 

Ion  and  electron  heating  associated  with  this  process  is  shown  in  Figure  5.3.  where  the 
ions  (left  panel)  and  electrons  (right  panel)  from  one  burst  are  shown  at  two  different  times 
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Figure  5.1.  Results  of  particle  simulations  of  two  simultaneous  expanding  plasma  clouds, 
showing  the  ions  of  each  burst  separately.  The  structure  is  more  pronoimced  on  the 
colliding  side.  (Lengths  are  normalized  to  c/u;,.) 
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Figure  5.3.  Results  of  the  same  run  showing  velocity  space  perpendicular  to  the  magnetic 
field.  Ions  are  energized  in  the  direction  of  the  other  cloud,  while  electron  heating  is 
isotropic.  (Velocities  are  normalized  to  the  expansion  velocity.) 
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{uLfit  =  4,  corresponding  to  the  appearance  of  structure  and  uJint  =  8  corresponding  to 
strong  penetration  of  the  two  debris  clouds).  There  is  strong  ion  heating  in  the  region 
of  the  overlap  of  the  two  clouds;  the  energetization  is  directed  toward  the  other  cloud. 
Electron  heating  is  also  quite  strong  and  essentially  isotropic. 

It  is  believed  that  the  heating  and  other  effects  associated  with  the  multiburst  case  are 
not  due  to  any  sort  of  ion/ion  streaming  instability.  Instead,  the  heating  is  attributed  to 
two  effects;  (1)  interaction  with  the  macroscopic  fields  and  (2)  interaction  with  the  waves 
generated  by  the  flute  instability.  The  first  process  involves  ions  from  one  burst  interacting 
with  the  macroscopic  electric  field  of  the  second  burst,  which  tends  to  accelerate  some  ions 
toward  the  center  of  the  second  cloud  (Figure  5.3).  There  is  also  the  effect  that  the  ions 
are  expanding  into  a  region  of  essentially  zero  magnetic  field  so  that  they  become  even  less 
“magnetized”  than  they  were  previously.  The  second  process,  namely  the  interaction  with 
waves  generated  by  the  instabihty  in  each  burst  separately,  leads  to  the  strong,  nearly 
isotropic  electron  heating.  Later,  in  Section  5.4  we  will  return  to  the  question  of  the 
possible  generation  of  direct  ion/ion  instabilities  zuid  their  possible  role  in  electron  amd  ion 
heating  in  multibursts. 

5.3  Hybrid  Simulation  of  a  Multiburst 

We  have  also  carried  out  a  few  multiburst  simulations  using  a  hybrid  (particle  ions,  massless 
fluid  electrons)  code.  The  calculations  are  similar  to  those  in  Section  1,  except  that  now 
there  are  two  equal  debris  clouds  expanding  perpendicular  to  a  background  magnetic  field. 
An  example  of  one  such  calculation  is  shown  in  Figmes  5.4-5. 7.  In  this  case  the  initial 
debris  density  of  the  clouds  relative  to  the  background  plasma  is  noo/no  =41,  and  the 
debris  mEisses  are  mo /mo  =  4.  The  separation  of  the  clouds  is  r,  =  Rg  =  12.8c/u>i.  The 
debris  expands  at  Mq  =  Fd/ua  =  0.5  with  a  directed  gyroradius  of  2c/u>,.  Figure  5.4 
displays  contour  plots  of  the  density  at  various  times.  One  sees  the  initial  expansion  of  the 
two  clouds  and  their  compression  into  thin  shells.  As  the  two  clouds  intersect,  there  is  a 
complex  collision  (fijf  =  32)  with  the  formation  of  one  cavity  by  =  48.  At  later  times 
the  clouds  recollapse  with  the  magnetic  field.  The  corresponding  contours  of  the  magnetic 
field  are  shown  in  Figure  5.5.  Again,  the  initizd  intersection  of  the  two  cavities  and  their 
eventual  merging  into  a  single  cavity  along  with  the  later  recovery  of  the  magnetic  field  is 
evident. 

A  time  series  showing  ion  dynamics  for  one  of  the  debris  clouds  is  presented  in  Figure  5.6. 
As  in  the  full  particle  runs,  one  again  sees  the  initi2il  radial  expansion,  the  plasma  being 
squeezed  into  a  thin  shell,  and  flute  modes  growing  on  the  surface.  (As  is  well  known,  in  the 
hybrid  simulations  of  VHANE  expansions  [Thomas  and  Brecht,  1988],  the  structure  forms 
somewhat  later  than  in  the  full  particle  simulations.)  As  in  the  single  burst  expeinsion 
(lower  right  psmel),  the  collapse  of  the  debris  carries  along  its  embedded  structure,  leading 
to  a  final  debris  cloud  with  r  ~  /2fl/2  that  remains  highly  turbulent.  However,  unlike  the 
full  particle  simulations,  the  debris  from  the  two  clouds  does  not  intersect.  In  this  case  as 
the  clouds  expamd,  air  ions  become  trapped  between  the  two  clouds,  giving  rise  to  a  large 
radial  electric  field  {Er  ~  Vpe)  that  tends  to  repel  the  debris  ions. 
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Hybrid  Simulations;  2-D  X  to  B 


Figure  5.4.  Hybrid  simulation  of  two  simultaneous  debris  expansions, 
contours  at  various  times. 


showing  ion  density 


magnetic  field 


(lit  =  0 


Qii  =  48 


0 


80  0 


80 


Figure  5.5.  Corresponding  magnetic  field  contours  for  the  same  case,  showing  the  merging 
of  the  two  clouds  to  form  one  large  cavity. 
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Figtirc  5.6.  Results  of  the  same  hybrid  run  showing  ions  &om  one  burst  in  the  x-y  plane 
at  various  times;  results  for  the  single  burst  case  at  the  same  final  time  shown  for 
comparison. 
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Corresponding  velocity  space  plots  of  these  same  ions  are  given  in  Figvire  5.7.  As  the  debris 
expands,  the  radial  velocity  decreases.  In  this  case,  however,  there  is  some  energization 
as  the  two  clouds  collide  {Qit  =  32).  The  more  energetic  ions  have  been  reflected  off  the 
interface  of  the  two  colliding  clouds.  As  in  the  runs  in  Section  1,  at  later  times  the  ions  axe 
heated  eis  the  magnetic  field  collapses.  Although  difficult  to  see  in  the  bottom  two  panels, 
the  ion  heating  in  the  multiburst  case  is  actually  greater  than  in  the  single  burst.  (While 
some  ions  in  the  single  burst  end  up  at  higher  velocities,  the  multiburst  case  has  a  flatter 
velocity  distribution  euid  hence  proportionally  more  ions  at  higher  energy.) 

5.4  Summary 

In  this  section  we  have  shown  examples  of  VHANE  multiburst  simulations  carried  out 
with  full  pairticle  and  hybrid  codes.  Both  types  of  calculations  exhibit  many  features  seen 
in  single  bursts  in  the  VHANE  regime,  especially  the  expansion  of  the  plasma  dragging 
the  magnetic  field  to  form  a  diamagnetic  cavity  and  the  development  of  structure  on  the 
surface  of  the  expanding  cloud.  And,  as  discussed  in  Section  1,  there  is  also  some  heating 
of  the  plasma  that  occurs  as  the  magnetic  field  recollapses.  However,  the  multiburst  cases 
also  show  new  and  interesting  effects  that  strongly  support  the  notion  that  multibursts 
cannot  simply  be  modeled  as  a  superposition  of  single  bursts.  These  include  additional 
heating  of  the  ions  and  electrons,  and  acceleration  or  deceleration  of  the  colliding  debris. 
A  first  analysis  of  these  simulations  suggest  that  the  new  effects  are  due  to  the  interaction 
of  the  plasma  in  one  burst  with  the  macroscopic  and  microscopic  fields  generated  by  the 
other  burst.  The  macroscopic  electric  field  that  is  produced  in  the  compression  shell  of  one 
burst  can  either  accelerate  ions  from  the  other  burst  toward  the  center  of  the  first  burst 
if  there  is  no  backgroimd  plasma,  or  reflect  the  ions  when  there  is  background  plasma 
compressed  between  the  two  clouds.  The  motion  of  ions  is  also  affected  by  the  absence  of 
the  magnetic  field  when  they  enter  the  cavity  of  the  other  burst. 

Microscopic  fields  due  to  the  lower  hybrid  drift  instability  that  generates  structure  on  the 
surface  of  the  clouds  also  contributes  to  the  ion  and  electron  heating.  In  particular,  it  is  felt 
that  the  large  and  isotropic  electron  heating  results  from  this  type  of  interaction,  rather 
than  from  the  excitation  of  some  sort  of  ion/ion  instability.  Such  cross-field  instabilities 
[Lampe  et  o/.,  1975]  can  occur  near  the  end  of  the  expansion,  when  the  relative  ion-ion 
speed  drops  below  the  local  Alfven  speed.  Such  instabilities  do  not  occur  in  the  hybrid 
formation  with  rrig  =  0  and  are  not  well  resolved  in  the  full  pairticle  simulations  (small 
milrric  precludes  mzmy  growth  limes  during  the  interaction).  Such  instabilities,  however, 
cannot  be  ruled  out  in  an  actual  VHANE  multiburst. 

VHANE  multibursts  would  be  a  very  interesting  study  using  the  AGNES  facility  at  NRL. 
The  VHANE  regime  is  easily  modeled  using  the  existing  experimental  conditions  [Ripin 
et  al,  1987].  The  above  simulations  suggest  that  there  are  several  new  effects  that  can 
be  investigated,  both  in  the  presence  or  in  the  absence  of  a  background  plasma.  In  ad¬ 
dition,  multiburst  experiments  could  be  used  to  help  resolve  the  old  issue  of  debris  air 
coupling,  by  examining  the  cross-field  ion  instabilities  mentioned  above  without  the  need 
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Figure  5.7.  Resxilts  for  the  same  nm  showing  the  ions  in  velocity  space. 


for  a  preionizer.  Laser  experiments  have  already  been  successfully  employed  to  study  ion- 
ion  interactions  for  the  ICF  program  [Berger  ei  al.,  1991]  and  similar  studies  could  be  very 
relevant  to  the  ezu-ly  time  problem.  We  are  presently  working  out  the  details  of  such  an 
experiment  and  will  present  the  results  in  a  future  report. 
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